UNCLASSIFIED 


AD  NUMBER 


AD810471 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  MAR  1967. 
Other  requests  shall  be  referred  to  Air 
Force  Weapons  Lab.,  Kirtland  AFB,  NM. 


AUTHORITY 


AFWL  ltr,  30  Nov  1971 


THIS  PAGE  IS  UNCLASSIFIED 


AFWL-TR-66-130 


AFWL-TR 

66-130 


CRITICAL  LOAD  CURVES  FOR 
REBOUND  CRACKING 

G.  S.  Samuelsen 
W.  Stuiver 
G.  R.  Abrahamson 

Stanford  Research  Institute 
Menlo  Park,  California 
Contract  AF  29{601}-6828 


/ 

TECHNICAL  REPORT  NO.  AFWL-TR-66-130 
March  1967 


AIR  FORCE  WEAPONS  LABORATORY 
Research  and  Technology  Division 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base 
New  Mexico 


AFVL-TF.-66-130 


Research  and  Tec'.mology  Division 
AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Conriand 
Kirtland  Air  Force  Base 
New  Mexico 


When  U.  S.  Government  drawings,  specifications,  or  other  aata  are  used  for 
any  purpose  other  than  a  definitely  related  Goverrroent.  procurement  operation, 
the  Government  thereby  incurs  no  responsibility  nor  ervy  obligation  whatsoever, 
and  the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any 
way  supplied  the  3aid  drawings,  specifications,  o-  other  data,  is  not  to  be 
regarded  by  implication  or  otherwise,  as  in  any  manner  licensing  the  holder 
or  any  other  person  or  corporation,  or  conveying  any  rights  or  permission  to 
manufacture,  uee ,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  thereto. 

This  report  is  made  available  for  study  with  the  understanding  that 
proprietary  interests  in  and  relating  thereto  uiil  not  be  impaired.  In 
case  of  apparent  conflict  or  any  other  questions  between  the  Government's 
rights  and  those  of  others,  notify  the  Judge  Advocate,  Air  Force  bysteas 
Command,  Andrews  Air  Force  Base,  Washington,  D.  C.  2C331. 

This  document  le  subject  to  special  export  controls  and  each  transmittal 
to  foreign  governments  or  foreign  nationals  may  be  made  only  with  prior 
approval  of  AF  (WLRP) ,  Kirtland  AF3 ,  N.F.  87117.  Distribution  cf  this 
document  is  limited  because  of  the  technology  discussed. 


DO  NOT  RETURN  THIS  COPY. 


RETAIN  OR  DESTROY. 


CRITICAL  LOAD  CURVES  FOR  REBOUND  CRACKING 


C.  S.  Samuel  sen 
W.  Stuiver 
C.  R.  Abrahjmson 


Stanford  Research  Institute 
Menlo  Park,  California 
Contract  AF29(601)-6828 


TECHNICAL  REPORT  NO.  AFWL-TR-66- 130 


This  document  is  subject 
to  special  export  controls  and 
each  transmittal  to  foreign 
governments  or  foreign 
nationals  may  be  made  only 
with  prior  approval  of  AFWL 
fWLRF) ,  Kirtland  AFB,  N.M. 
Distribution  of  this  document 
is  limited  ueuauae  t  f  the 
technology  discussed. 


An.^-TR-66-130 


FORFWnRD 


This  report  was  prepared  by  the  Pouller  Laboratories,  Stanford  Research 
Institute,  Menlo  Park,  California  under  Contract  AF29(601)-6828 .  The  research 
was  performed  under  Program  Element  6.16.46.01D,  Project  5710,  Subtask  15.014, 
ar.d  was  funded  by  the  Defense  Atomic  Support  Agency  (DASA)  . 

Inclusive  dates  of  research  were  June  1965  through  September  1966.  The 
report  was  submitted  13  February  19b7  by  the  AFVL  Project  Officer,  Lt  Walter  D. 
Dittir.er  (WLRP)  .  Contractor's  report  number  is  FGU-5569. 

Technical  direction  was  provided  by  Lt  R.  C.  Brightman  and  Lt  W.  Dittmer  in 
coordination  with  Lt  L.  T.  Montuili,  all  of  the  Physics  Branch  of  the  Air  Force 
Weapons  Laboratory  (WLRP) . 

The  authors  are  indebted  to  Betty  Bain,  who  provided  invaluable  assistance 
with  the  numerical  calculations  presented  in  Section  2,  and  to  Leonard  Dary,  who 
was  responsible  for  the  experimental  work. 

This  report  has  been  reviewed  and  is  approved. 

WALTER  D.  DITTMER 
lit.  l!SAF 
Project  Officer 


EDGAR  M.  MbNYON 
iJSAF 

Chief,  Physics  Branch 


Colonel,  USAF 

Chief,  Research  Division 


ABSTRACT 


It  is  shown  that,  under  an  external  load,  a  cylindrical  pl»stlc 
shell  mounted  on  a  steel  mandrel  can  rebound  ut  a  ve.ocity  sufilcient 
to  cause  cracking  of  the  shield  due  to  hrop  tension.  Using  numerical 
alculations,  critical  load  curves  are  found  for  the  combinations  of 
peck  pressure  and  impulse  which  divide  the  loads  that  produce  cracking 
from  those  that  do  not.  These  curves  are  found  to  be  consistent  with 
semigraphlcal  solutions  based  on  the  method  of  characteristics.  The 
latter  method  is  also  used  to  obtain  critical  load  curves  for  spall. 

Experiments  that  verify  the  theory  of  rebound  cracking  are 
reported. 
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SECTION  1 


INTRODUCTION 

As  illustrated  in  Fig.  1.1,  a  cylindrical  plastic  Shield  supported 
by  a  metal  base  structure  may  fail  by  rebound  cracking  under  a  suddenly 
applied  surface  load.  In  this  report  the  shield  is  taken  to  be  thir 
compared  to  its  radius  and  compared  to  the  thicknear  of  the  supporting 
base  structure;  the  load  is  taken  as  uniform  over  the  surface.  Under 
these  conditions,  the  rebound  velocity  of  the  shield  can  be  found  using 
one-dimensional  wave  theory.  To  calculate  the  hcop  st.  *as  developed, 
the  subsequent  motion  of  the  shield  is  taken  to  he  a  linear  oscillation. 
Cracking  occurs  if  the  hoop  stress  developed  during  the  outv.ard  motion 
exceeds  the  strength  of  the  shield. 

For  loade  with  a  step  rise  and  gradual  decay,  the  combinations  of 
peak  pleasure  and  total  impulse  which  result  in  cracking  would  be  ex¬ 
pected  to  appear  as  shown  in  Fig.  1.2.  For  loads  of  duration  less  than 
twice  the  transit  time  through  the  shield,  the  outward  displacement 
depends  only  on  me  impulse.  As  load  duration  increases,  the  peak  pres¬ 
sure  must  be  increased  to  maintain  the  rebound  velocity.  Moreover,  for 
the  loads  acting  after  separation,  a  higher  rebound  velocity  is  necessary 
to  reach  the  critical  displacement.  Hence,  as  load  duration  increases, 
peak  pressure  (and  uius  impulse)  required  to  cause  cracxing  increase  as 
Indicated  in  Fig.  1.2. 

In  the  present  investigation  methods  for  determining  cracking 
boundaries  Are  developed  and  the  theory  of  rebound  cracking  is  checked 
experimentally.  Cracking  boundaries  are  determined  for  a  typical  plastic 
suppcnc.1  by  a  base  of  thick  aluminum,  a  base  of  thin  aluminum,  and  a 
base  of  thin  aluminum  over  a  rigid  foam  for  pulses  with  a  step  rise  and 
exponential  decay.  Moie  limited  calculations  are  made  for  pulses  with 
a  step  rise  and  linear  decay.  Rebound  velocities  are  obtained  from  nu¬ 
merical  calculations  based  on  the  -thod  and  the  maximum  displacement 
is  found,  regarding  the  plastic  cylinder  as  a  linear  oscillator.  The 


1 


FiG.  1.1  EXPERIMENT  DEMONSTRATING  REBOUND  CRACKING.  The  explosive  drives 

.U:#,U  -  -♦  «i_  .»..i  c  ...1  t  I  I  I  I 

-  .  v  ^  uyw..,  L-i  me  mcci,  omvc  i * i tr  bieeri  nu:  a  mgner  shock  impedonce, 

tne  sh'eld  rebounds.  If  the  rebound  velocity  is  sufficient,  crocking  occurs  os  rhcwn. 
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FIG.  1.2  SHAPE  OF  CRACKING  BOUNDARY 


cracking  boundary  in  the  pressure-impulse  plane  is  the  line  separating 
the  region  of  cracking  from  the  region  of  no  cracking.  The  cracking 
boundary  is  also  determined  using  a  graphical  method  based  on  lirear 
v v c  theory .  Thu  boundaries  obtained  the  t— ^  meuiuua  aic  iuunu  i,o 
be  consistent. 

The  graphical  technique  is  also  used  to  obtain  the  spall  boundary 
in  the  pressure-impulse  plane.  The  boundary  obtained  is  consistent  with 
the  results  of  the  numerical  calculations. 

A  limited  series  of  experiments  were  undertaken  to  check  the  valid¬ 
ity  of  the  theory  of  rebound  cracking.  The  experimental  results  agree 
with  the  predictions  within  experimental  error. 

In  Section  2  tne  numerical  calculations  are  described  end  the 
cracking  boundaries  obtaired  are  plotted.  In  Section  3  an  introduction 
to  the  graphical  analysis  is  given  and  the  results  arc  compared  with 
those  of  the  numerical  calculations.  A  detailed  description  of  the 
graphical  method  is  given  in  Section  4.  Tha  experiments  undertaken  to 
verify  the  theory  are  described  in  Section  5. 
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SECTION  2 


NUMERICAL  SOLUTIONS 


2. 1  Introduction 

The  response  of  a  shield  on  a  thick  base  has  been  extensively 
analyzed  for  uniform  [  1  ]  and  exponential  pulses  (Section  4)  for  lir.ear- 
clastic  materials,  using  the  method  of  characteristics.  However,  under 
high  pressures  some  materials  exhibit  strong  nonlinear  properties  and 
respond  plastically.  Such  conditions  introduce  complex  interactions 
between  variables  ar.d  necessitate  a  different  theoretical  approach. 

In  the  present  study  a  computer  code  was  modified  to  study  the 
effects  of  pulses  with  a  step  rise  and  exponential  or  linear  decay, 
applied  uniformly  to  the  exterior  of  a  cylindrical  model  for  which  the 
modes  of  failure  are  cracking  and  spalling  of  the  shield.  Three  models 
were  investigated:  (1)  two -layered  composite  with  a  thick  base,  (2)  two- 
lavered  composite  with  a  thtn  base,  and  (1)  t hree— 1 avered  composite  with 
a  thin  base  backed  by  foam. 

The  primary  purpose  of  this  work  was  to  determine  critical  load 
boundaries  for  rebound  cracking  of  the  shield,  taking  account  of  non¬ 
linear  and  plastic  effects  on  stress-wave  propagation.  The  calculations 
are  divided  into  two  parts:  (!)  early-time  stress-wave  propagation  in 
the  composite  prior  to  separation,  and  (2)  late-time  structural  dis¬ 
placement  of  the  shield  after  separation.  The  stress-wave  calculations 
yield  the  rebound  velocity  of  the  shield,  which  is  the  initial  condition 
for  the  structural  displacement  calculations. 

The  shield  is  taken  to  be  thin  relative  to  the  radius,  permitting 
a  one-dimensional  model  to  be  used  for  the  stress-wave  calculations. 

The  stress-wave  calculations  are  made  using  the  Q-method  of  von  NeumAnn 
e.nd  Richtmyer  [2],  which  has  proven  very  useful  for  nonunlform  wave 


flow,  particularly  where  nonlinear  material  properties  are  involved. 

The  method  Is  based  upon  finite  difference  equations  and  an  artificial 
viscosity  function  which  removes  discontinuities  in  the  flow  field; 
wave  interactions  are  handled  autoniati  ally  and  it  is  unnecessaiy  to 
locate  the  wave  fronts  as  with  the  method  of  characteristics.  The  code 
used  here  is  based  on  a  code  developed  by  White  and  Griffis  l  3 J .  Wood  [ 4 J , 
and  Morland  [5],  and  generalized  by  Lrkman  l 6 ] . 

The  late-time  structural  model  used  is  a  ring  in  uniform  radir.l 
motion,  regarded  as  a  linear  oscillator.  The  maximum  displacement  is 
found  for  the  rebound  (initial)  velocity  and  the  incident  load  acting 
after  separation.  Cracking  is  taken  to  occur  11  the  displacement 
exceeds  a  critical  value. 

""he  three  models  investigated  are  described  in  Table  2.1.  For 
all  models,  the  shield  and  base  were  epoxy  and  aluminum  and  the  shield 
thickness  was  2.54  cm.  Model  1,  the  two  layered  (thick  base)  composite, 
was  the  one  most  extensively  analyzed.  For  this  model,  the  influence  of 
linearity  and  plasticity  were  studied  for  exponential  and  triangular 

,  —  —  -  r  v  1 - 

pu  1  O  V,  pACyXXAVO* 

Models  2  and  3  were  used  to  study  the  effect  of  rigid  foam  as 

a  third  layer.  The  thickness  of  the  thin  base  va:  0.635  cm.  For  Model 

.  3 

3,  a  thick  layer  of  7  lb/ ft  foam  was  added  to  Model  2. 


Table  3.1 


kATTJUA L  PROPISTIES 


Material 

Dan  si  t  > 

(g  c»3 ) 

Sound  Spaad 
(c*/>a«c  ) 

Acouitlc 

I spadarica 
( btxi«c/ca) 

Yl  aid 

Strap* 

Ubmr) 

Pol saon ' • 
Ratio 

Modal 

Thioknaaa 

(Cfc) 

Kean 

Radiua-r 
(cn)  * 

Epoxy  C-7 

118 

0.265 

0.313 

3.0 

0.  4 

2 . 54 

36.  1 

A  lust  rum 

2.  «fi 

0.640 

1.  TS1 

5.0 

0.3 

1 

inf  tm  ta 

18081-T8) 

2,3 

0.635 

Rigid  foaa 

O.il 

(t  lb/ it  ) 

0.253 

0.028 

_ 

o.on 

3 

”inf lnl te 

la  Section  2.2  the  theoretical  bases  for  the  stress-wave  and 
Structural  displacement  calculations  arc-  given  and  the  failure  criteria 
are  described.  In  Section  2.3  the  critical  cracking  boundaries  deter¬ 
mined  from  the  numerical  calculations  are  presented. 

2 . 2  Approach 

2.2.1  Basic  Equations  for  Stress-Wave  Calculations 


The  shield  and  base  materials  arc  represented  by  the 
e lust op  la  si ic  model  ot  Fig.  2.1.  fcith  this  model,  the  material  behaves 
elastically  along  ba,  in  which  a  is  the  yield  point  in  compression. 
Segment  ae  is  the  locus  of  plastic  loading  states,  ef  the  locus  of 
elastic  unloading  stages,  and  fb  the  locus  of  plastic  unloading  states. 
The  mean  pressure  and  maximum  resolved  shear  stress  are  defined  by 

dp  =  |  dfp^  +  py  -  p^)  (2.1) 

up  “  UU 

di  =  -■*  2  ~-V  (2.2) 

and  the  compressive  stress  in  the  direction  of  propagation  is  given  by 

dp  =  dp  -i  ^  d  r  (2.3) 

x  3 

The  relationship  between  the  mean  pressure  and  density  (Hugonict)  was 
assumed  to  be  of  the  form 

p  =  AH  +  (2.4) 

where  n  =  p/ o  -  1  and  A,  B,  and  C  are  constants  depending  on  the 
material  (see  Fig.  2.1). 

Equation  (2.3)  is  the  basic  relation  used  in  the  numerical 
calculations  and  was  specialized  for  the  various  segments  of  Fig.  2.1. 
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FIG-  2.1  ELASTOPLASTIC  CYCLE 


For  the  elastic  region,  segment  ab  ,  (2.3)  becomes 


Up  =  (K,  +  r  G)  dr', 

x  b  j  x 


(2.5) 


where  Is  toe  bulk  modulus  and  G  the  rigidity  modulus.  Using 


K  -  £  (._!__)  =  .dP 
b  3  V1  -  2V  dV 


(2.5a) 


G  = 


2(1  -  'J ) 


(2.5b) 


dc  = 

C 


(2.5c) 


(2,5)  can  be  written  as 


Taking  v  9 s  constant ,  (2.6)  Integrated  to 


_  ,  1  -  v  - 

p  -3  —  p 

xc  1  -  v 


(2.7) 


»hci'i  the  subscript  r  depot r*  the  clastic  segment  nb  .  Foi  uniaxial 
strain,  the  compressive  stress  normal  to  the  direction  of  wave  propaga¬ 
tion  is  given  bv 


dp 


(K  -  2p/ 3 )dt 


(2.8) 


and  using  (2.2)  and  (2.6),  the  resolved  slieai  strers  becomes 


3  -  ,1  -  2v. 

T  =  2  p  <TTV) 


(2.9) 


For  the  Treses  and  the  von  Ml  sea  yield  conditions  in  this  geometry,  the 
maximum  resolved  shear  stress  is  determined  by  the  yield  stress  Y  =  2t. 
Consequent ly . 


Y  =  3p  (,-  ;  ?>) 


(2  10) 


Along  ae  ,  (2.3)  Integrates  to 


P  *  P 
Xl  xa 


P  -  P  —  ( Y  -  Y  ) 
a  3  a 


(2.11) 


where  the  subscript  t  denotes  loading.  Along  fb  , 


xut 


■  pvf  =  p  ■  p,.  -  r  ;y  -  v  > 
xl  1  o  i 


(2.12) 


where  u-1  uviivtr-.  unloading.  me  other  segments  are  determined  in  a 
similar  manner. 


Foi  the  foam  (Model  3),  the  dependence  of  pressure  upon 
specific  volume  Is  shown  In  Fig.  2.2.  For  the  elastic  states,  the  com¬ 
pressive  stress  in  the  direction  of  propagation  is  given  by 


6 


WfSSv^*'. 


As  the  material  is  compressed  below  V ^  ,  the  pressure  is  given  by 


a'  (V .  -  V  , 


(  -  MVj  -  V)  +  v^]2 


(2.15) 


The  equations  of  continuity  and  motion  for  cne-dimensional 
fi  ns  In  Lagrangian  coordinates  nre 


dV 


5u 


co  dx 


(2.16) 


cu 

At 


d(p  +•  4') 

x 

dx 


(2.17) 


where  u  is  particle  velocity,  and  Q  is  the  artificial  viscosity  param¬ 
eter.  In  the  code  used  here,  Q  was  taken  to  be 


■  v  k  (!  I*  I  ‘2  *  “■ ) 


(2.181 


where  i  =  1,7,  ti  =  0.2,  and  c  is  the  local  sound  speed. 

For  the  calculations,  (2.16)  sad  (  ’.17)  were  written  in 
difference  form,  programmed  in  ALGOL,  ana  integrated  for  the  Alven 
initial  conditions  using  the  relations  of  Figs.  2.1  and  2.2  for  the 
elastoplastlc  ar.d  foati  layers.  The  program,  is  characterized  ay  the 
following  progression  of  3teps: 


B. 


For  each  time  step,  the  particle  velocity  is  found 
in  every  cell ; 

the  di splaceme-nt  of  ce.ch  cell  is  determined  and  the 
corresponding  value  of  P-  and  V, 


for  me  shield  and  base,  p  i»  determined  as 

X 

follows: 


1,  If  p  is  greater  than  that  of  the  previous 
time  ctep,  loading  occurs  and 

a.  if  the  pressure  p  fron  (2.7)  is  nbeve 
p  ,  from  (2.11)  tft£  value  of  p  is 
l  aiten  ;  otherwios 
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1 


\ 


I 


b.  the  pressure  i.6  la  the  elastic  region  and  the 

value  of  p  is  taken, 
xe 

2.  It  P  is  less  than  that  o i  the  previous  time 
step,  unloading  occurs  and 

a.  if  pxfc  is  below  Pxu£  from  (2.12)  the 
pressure  is  equated  to  Pxu^  .  otherwise 

b.  the  pressure  is  in  the  elastic  unloading 
region  end  the  value  of  pxe  is  taken. 

D.  For  the  foac ,  if  included,  px  is  determined  from 
(2. 13 )through(2. 15) .depending  upon  the  value  of  the 
specific  volume  V  . 

The  output  of  the  calculations  includes  (1)  the  occurrence 
and  tifc*e  of  spall,  and  (2)  the  time  of  separation  ana  rebound  velocity 
of  tne  shield.  The  latter  becomes  the  initial  velocity  for  the 
structural  displacement  calculations. 


2,2.2  Structural  Displacement  of  the  Shield 

As  outlined  above,  the  structural  model  is  taken  as  e 
j  ir.p  behaving  as  »  linear  oscillator.  The  equation  of  motion  for  e 
linear  oscillator  (neglecting  damping)  is 


k 

+  —  x 

Si  El 


F(t ) 


(2.19) 


where  sn  Is  the  mass,  F(t)  the  forcing,  function,  and  k  the  spring 
constant.  la  the  preaent  problem,  h/u  la  related  to  tne  hoop  modulus 
Eg  by 


v 

IE 


cr 


m 


(2.20) 


in  which  r  la  the  mean  radius  of  the  shield.  Since. 


(?.21) 


-. — 'VSku't*^.  j*;aaoiMuS.i 


12 


and 


c 

r 


(2.20)  may  be  written 


q 


m 


*- 


y 

2v 


1 


Using  (2.20)  and  (2.23),  the  solution  to  (2,19)  may  be  written 


x 


or 


x 


where 


and 


x  cos  q  t  +  —  sin  q  t 
s  m  t  m 

m 


1_ 

q_ 


J  F(t)  sin  q^(t  -  t  x)dt  ' 


v 

x  cos  q  i  •»  —  sin  o  t  +  G  cos  q  t  +  H  sin  q  t 
a  m  q  m  m  n. 

a 


G 


1 


+ 

\  F(t )  sin  q  t  'dt 
J  m 

o 


H 


F(t )  cos  q  t  'at  ' 
n* 


x  =  initial  displacement 

s 

t  =  time  of  separation 

n.  -  mass  of  shield  =  prr( x ^  -  r^)  (pe~  unit  length) 

r  -  outer  redius 

o 

r}  =  inner  radius 

F(t)  -  P(t)  x  area  =  P(t)  2n(r  +  —  +  x)  (per  unit  length) 

m  2 


(2.22) 


(2.23) 


(2.24) 


(2,24a) 


(2.24b) 


(2.24c) 
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For  exponential  and  triangular  pulses  respectively,  we 


have 


Pvt)  «■  PQe  At  (Fig.  2.3) 


P(t)  =  P  (1  -  h  (Fig.  2.4) 
O  T 


FIG.  2.3  EXPONENTIAL  PULSE  FIG.  2.4  TRIANGULAR  PULSE 


Substituting  these  into  the  above  equations  end  Integrating  yields 


x  = 


x  cos  q  t  +  —  filn  q  t  +  5F 
a  m  q  ta 

m 

Y~z-gf 


(2.25) 


where 


P  B  (r  ♦  £> 
o _ at  2 

pch 


S 


f  2, 25a) 


and 


x2  +  a 


2  i  V  -  " 


-\t 


m 


-  X  sin  q  (t  -  t  )  +  q  cos  q  (t  -  t  ) 
s  m  m  s 


C  exponential ) 


(2.250) 


or 


F  =  — 
q_. 


1  -  |  -  {=S  -  1) 

T  T 


cos  q  (t  -  t  )  +  — -  sin  q  (t  -  t  ) 
m  s  q^T  m  s 


(triangular) 


(2.25c) 


2.2.3  Failure  Criteria 

Spall  Stress.  Spall  failure  was  assumed  to  occur  in  the 
epoxy  shield  if  the  tensile  stress  exceeded  the  yield  stress  of  the 
epoxy,  3. 0  kbar . 

critical  Velocity  for  Cracking — Short  Pulses.  If  the  inci¬ 
dent  pul3e  does  not  act  beyond  separation,  a  critical  separation  velocity 
suificient  to  develop  hoop  stresses  to  the  point  of  cracking  may  be  de¬ 
fined.  The  reaicual  velocity  at  separation  is  the  average  particle 
velocity  across  the  layer  at  time  of  separation.  If  vg  is  the  resi¬ 
dual  velocity  at  separation,  the  kinetic  energy  (per  unit  area)  is 


KE 


1  .  .  ,2 

2  ph  'V 


(2.26) 


For  uniform  radial  displacement,  the  strain  energy  is 


iLc2 

2E«  e 


where  E  is  the  modulus  in  the  circumferential  direction.  Equating 
the  strain  energy  (2.27)  to  the  kinetic  energy  at  separation  yields 

E . 


For  the  critical  condition, 
Using  (2.26)  and  solving  or 


a  becomes  the  maaimur.  hoop  stress  n 

6  0cr 

vg  yields  the  critical  separation  velocity 


v 

cr 


(2.29) 


Using  (2.21)  and  (2.22),  (2.29)  may  be  written 


cr  pc 


(2.30) 


Taking  oQ  ■  3  kbar, 
ocr 


v  =  12.88  cm/ ms 

cr 


Critical  Displacement  for  Cracking — Long  Pulses.  Pulses 
which  act  after  separation  must  be  accounted  for  in  calculating  the 
structural  displacement  of  the  shield.  From  Hooke's  law 


r  -  r 

Si 


r 

m 


°0E0 


(2.31) 


where  r  is  the  radial  displacement  of  the  mean  radius  r  .  Using 

m 

(2.21)  and  (2.22),  and  putting  £x  ■  (r  -  r  ),  a  =  <y  ,  (2.31) 

cr  cr  m  bcr 

yields 


Ax  *>  ^ — -  r 

cr  ©cr  oc  m 

r 


(2.32) 


For  the  properties  of  Table  2.1  the  critical  displacement  for  the  epoxy 
shield  Is 


Axcr  •  2.48  cm 
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2.3  Results 


2.3.1  Introduction 

The  results  are  presented  in  two  parts.  First  the  results 
for  the  two-layered  composite  with  a  thick  base  ere  given;  these  show  the 
effects  of  differences  in  the  constitutive  relations  for  an  exponential 
pulse  and  comparing  the  effects  of  exponential  and  triangular  pulses. 
Next,  the  results  are  given  for  the  two-layered  composite  with  a  thin 
base  and  the  three-layered  composite  with  a  thin  base  and  a  thick  foam 
layer,  all  for  an  exponential  pulse. 

2.3.2  Moccl  1 — Two-Layered  Composite  with  a  Thick  Base 

Exponential  Pulse.  Critical  load  boundaries  for  rebound 
cracking  of  the  shield  were  generated  for  an  exponential  pulse  (Fig.  2.3) 
applied  to  an  epoxy  shield  and  aluminum  base  described  by  the  model  of 
Fig.  2.1.  Effects  of  differences  in  material  properties  were  examined 
for  the  following  cases: 

1.  Nonlinear-Nonelastic  [a,  B,  C  /  0  in  (2.4)] 

2.  Nonlinear-Elastic  (p  =  in  Fig.  2.1) 

xa 

3.  Linear-Elastic  (B,  C  =  0;  p  =  «°) 

4.  Linear-Nonelastic  (B,  C  =  0) 

The  differences  between  the  nonlinear  and  linear  constitutive  relations 
can  be  seen  in  Fig.  2.0. 

Figures  2.6a  end  2.6b  show  the  results  of  the  computer 

* 

runs  and  the  cracking  boundaries  determined  from  them.  As  expected, 
the  boundaries  all  exhibit  a  minimum,  are  impulse-dependent  for  low 
impulses,  and  rise  monotonically  beyond  the  minimum.  Spall  occurred 
as  shown. 

* 

See  Tables  2.2  through  2.5  for  numerical  values. 
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FIG.  2.5  COMPARISON  OF  NONLINEAR  AND  LINEAR  EQUATIONS 
OF  STATE  FOR  EPOXY 


Effect*  of  Differences  In  Constitutive  Relations  of 
Shield  and  Base.  Eff— of  differences  In  the  constitu¬ 
tive  relations  of  the  shield  and  base  are  shown  In  Fif .  3.6c  (obtained 
from  2.6a  and  b).  For  the  range  of  differences  considered  here,  non¬ 
linear  effects  (l.e.,  B,  C  t  0)  are  much  more  significant  than  plasticit 
effects. 


IB 


FIG.  2.6  (Concluded)  CRACKING  BOUNDARIES  FOR  EXPONENTIAL  PULSES 
(c)  (o)  end  \b)  combined 

Nonlinear  eifects  result  In  a  dilference  c f  about  40£  In 
critical  pressure  at  the  minimum.  To  the  right  of  the  minimum,  nonlinear 
effects  produce  differences  of  about  40%  or  leas  in  both  pressure  and 
impulse.  To  the  left  of  the  minimum  tin  the  region  where  the  curves 
rise  rapidly)  only  the  Impulse  is  significant  and  differences  in  Impulse 
due  to  nonlinear  effects  are  about  15%.  In  contrast  to  the  nonlinear 
effects,  plasticity  effects  result  in  differences  of  less  than  10%  every- 
%utr6 , 


Nonlinearities  affect  the  critical  load  curves  prima:.-ily 
through  increases  in  shock  velocity  with  pressure,  which  result  in 
earlier  separation  for  the  nonlinear  material.  This  decreases  the  re¬ 
bound  velocity  and,  due  to  the  higher  load  remaining  st  the  earlier 
separation,  causes  the  shield  to  experience  more  rapid  deceleration 
after  separation.  Thus,  the  load  amplitude  required  to  develop  a  given 
displacement  will  be  higher  for  the  nonlinear  material. 


30 


T»bl»  7.2 

EXFOfllUrriAL  PUUiX  <HOdU*lAA-HO*rt LACTIC  ) 


F 

1 

<#»••€  ~l) 

I 

(10®  t«Di) 

Rue 

T1m  of 
Saparat loe 
(Faae ) 

— 

l«ptrn  ion 
V«1«11  r 

(Cf'  MK  > 

— 

lUllMUK 

Dlftplacaaent 

(cal 

Crack 

Spall 

Tlae 

(naec  ) 

— 

13.2 

0.  076 

22 

■ 

_ 

14.9 

■ 

mm 

13  28 

0.4 

2.00 

93 

»fii.  J 

- 

19.9 

- 

non# 

0.33 

3. 00 

P4 

20 . 00 

- 

11.4 

- 

■ 

DO 

non# 

4  8 

0.  083 

90 

20.38 

- 

is. ; 

»#• 

12  39 

c 

•a 

c 

0.  308 

1.30 

69 

13.82 

- 

38.8 

>•• 

13  74 

0.K87 

1 . 60 

91 

13.67 

- 

21.2 

- 

2.33 

>*■ 

nona 

BSeM 

2.30 

0.08 

31 

31.01 

- 

13.  n 

- 

2.31 

>#• 

14.  18 

isRSylM 

2.3u 

0.  087 

49 

21.33 

- 

14.0 

2.70 

>*• 

13.96 

■  1 

o.2* 

0.87 

4a 

14.67 

- 

37.0 

4.94 

nona 

hzb 

0.17 

1.176 

30 

20.  19 

- 

3.0 

- 

1.31 

no 

nana 

0.06 

1.0 

0.  06 

42 

23.24 

- 

11.0 

2.  12 

no 

16. 16 

0  06 

0.80 

0.  075 

4^ 

22  88 

- 

13.0 

2. 32 

>*■ 

16.00 

0.06 

0.  70 

0.  086 

43 

22.  28 

- 

14.6 

2.63 

>•• 

13.63 

o.(y\ 

0.13 

0.  40 

43 

17.  33 

- 

3.67 

>•• 

nona 

0.08 

0.  13/ 

0.438 

87 

17.47 

- 

3.30 

yea 

nona 

0.06 

0.  13 

0.-261 

46 

17.90 

- 

1.83 

no 

nona 

0.06 

0  12 

0.30 

** 

18.  71 

- 

0.80 

no 

nona 

0.03 

0.  40 

0.073 

36 

23.30 

- 

2.64 

ye» 

17.51 

0.03 

0  IS 

0.20 

38 

20.83 

- 

3  67 

ye# 

nona 

0.  03 

0.  12 

0.23 

40 

2J.47 

- 

21.3 

2.81 

y». 

nona 

0.03 

0.  10 

0.30 

37 

20.24 

■ 

BSB 

1 .  Jo 

no 

nona 

0.013 

0.325 

0.046 

3* 

25  31 

■ 

wxM 

1  87 

no 

20.  98 

0.  013 

0.20 

0.073 

35 

24  20 

■ 

sis 

2.35 

dj» 

12.  53 

0.015 

0.  15 

0.  10 

33 

S3  83 

14.0 

- 

2.52 

nona 

0.013 

J.  10 

0.13 

32 

33.15 

14.0 

- 

1.71 

r-enc 

- - 

_ 

_ 

. 

- L 

T»bl»  2.3 

EXPOWWTIAL  ?IUE  (UKXAR-EIASTIC  ) 


P 

<m>ar) 

\ 

(Haac  ) 

I 

<10€  taps) 

Bun 

Tiaa  of 
Separation 
(^•ec ) 

Saparat  Ion 
Valjyl  tv 
<  c*/»aac  ) 

M**li*ia 

Di  aplocaaant 
(ca) 

spii: 

Tip* 

(MB#C  ) 

mm 

3.  3 

0.  037 

" 

55 

19.00 

- 

11.0 

• 

2.12 

no 

in 

3.0 

0.067 

37 

26.45 

- 

13.6 

- 

2.40 

cloaa 

1&1 

0.  137 

1.66 

83 

19. Ov 

- 

68.0 

- 

3.16 

y#a 

m 

0.20 

0.  13 

1.54 

62 

19.  78 

- 

52.3 

- 

2 . 07 

no 

ifci 

0.20 

0.  n 

1*2 

38 

23  ii 

- 

13.1 

- 

0.  oo 

no 

nana 

0.  Oj 

0.  30 

0.06 

33 

27  94 

- 

12.2 

- 

3.30 

no 

16.60 

0.  03 

0.073 

0.  40 

34 

26.02 

- 

39.0 

- 

3 . 62 

23.95 

0.03 

C.G6 

0.30 

36 

23.78 

- 

38.0 

- 

1.17 

nona 

0.013 

0.22 

0.068 

60 

2d.  85 

- 

13.7 

- 

2.64 

20.93 

0.013 

0.07 

55 

28.42 

- 

23.0 

- 

2.48 

.6.33 

0.01 

0.  10 

32 

29.44 

- 

15.2 

- 

2.  37 

24.57 

0.01 

0  073 

wmm 

29.39 

“ 

16.3 

~ 

2.  13 

wBM 

28. u4 

Triangular  Pulse.  :  he  cracking  boundary  for  a  triangular 
pulse  for  the  nonlinear,  nonelastic  case*  Is  given  in  Fig.  2.7  together 
with  that  for  an  exponential  pulse.  The  curves  coincide  in  the  impulse- 
ssnsitive  region,  as  expected.  The  minimum  for  the  triangular  pulse  is 
about  50%  below  that  for  the  exponential  pulse.  To  the  right  oi  the 
minimum  the  curves  difier  by  about  30%. 


F'G.  2.7  CRACKING  BOUNDARIES  f  OR  EXPONENTIAL  AND  TRIANGULAR 
PULSES  TOR  NONLINEAR.  MONFLASTIC  CASE 


Ste  Table  2.6  for 


numeri co 1 


values , 
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The  triangular  pulse  differs  from  the  exponent i al  in  that, 
for  the  same  pent  pressure  and  total  impulse,  the-  exponential  pulse  has 
a  higher  initial  decay  rate  and  a  longer  duration.  Thus,  for  triangular 
and  exponential  pulses  of  the  satr.e  pea':  pressure  in  the  region  where 
curves  rise  to  the  right.,  less  reflected  momentum  (Irons  the  base)  is 
realized  for  the  expor  -n^ial  pulse  and  more  of  the  pulse  acts  after 
reparation.  Thus,  to  uevelop  the  same  maximum  displacement  for  a  given 
peaK  pressure,  the  exponential  pulse  requires  a  shorter  duration  and 
hence  less  momentum. 

Above  0.05  Mbar  in  the  region  where  the  curve  rises  to  the 
right,  the  numerical  results  suggest  that  0  distinct  break  exists  in 
the  critical  boundary.  This  le  different  from  the  exponential  results 
which  show  a  eontir.UG'dSl.v  ilbing  curve.  The  dashed  line  is  the  ivur.d&*y 
a:  it  would  be  expected  to  develop  for  Increasing  values  oi  impulse. 

2.3.3  Model  2 — Two-Layered  Composite  with  a  Thin  Base 

To  provide  a  basis  for  evaluating  the  effect  of  a  third 
layer  oi  roam,  s  few  cases  were  run  without  the  foam.  The  mode)  is 
described  in  Table  2.1.  The  calculations,  employing  the  nonll near¬ 
nonelastic  constitutive  relation  and  exponential  pulse.  Indicated  that 
the-  rarefaction  from  the  free  surface  of  the  base  caused  greater  Inward 
particle  velocities  to  be  realized  in  tne  base  than  in  the  shield.  As 
a  result,  the  base  separated  from  the  shield  when  the  rarefaction  arrived 
at  the  intei  Iftce.  Because  the  average  velocity  of  the  base  was  over  three 
times  larger  than  that  of  the  shield,  it  w.-.s  assumed  that  the  shield  would 
not  re-impact  the  base  during  the  inward  notion.  Alter  separation,  the 
snielc  motion  was  regarded  as  ■-  linear  oscillator.  As  before,  cracking 
was  taken  to  occur  if  the  outward  displacement  exceeded  2.4S  cm.  Al¬ 
though  the  Inward  displacement  was  usually  equal  to  or  greater  than  the 
outward  displacement,  failure  in  tension  was  assumed.  This  is  reasonable 
for  a  material  with  a  yield  strength  in  compression  which  exceeds  that 
in  tension 
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The  results  are  Indicated  111  Fig.  2.7  by  the  points 
labeled  322,  323,  and  324.*  Although  all  three  points  are  within  the 
cracking  region  for  a  thick  base,  only  324  exhibited  cracking. 

2,3.4  Model  3-~Three-Layered  Composite  (Thin  Base  with  Foam) 

For  thin  model  the  thin  base  used  for  the  calculations 
described  immediately  above  was  backed  with  a  thick  layer  of  rigid  poly- 

3 

urethane  loom  of  density  7  lb/ft  (see  Table  2.1  and  Fig.  2.2).  The 
foam  would  be  expected  to  affect  the  results  by  attenuating  the  rare¬ 
faction  wa va  iron  the  back  surface  of  the  base. 

The  results  of  the  calculations  are  Indicated  in  Fig.  2.7 
by  the  point  labeled  304.  The  maximum  outward  displacement  was  2.25  cm. 
For  the  same  load  (point  322),  the  thin  base  without  foam  gave  a  maximum 
outward  displacement  of  2.2S  cm.  Thus,  the  foam  did  net  significantly 
affect  the  response. 

A  foam  of  greater  density  might  have  prevented  the  base 
iron  separating  as  a  result  of  the  free  surface  rarefaction.  This  would 
have  increased  the  time  during  which  the  foam  could  influence  the  mo¬ 
mentum  exchange  with  the  shield  and  the  affect  on  the  velocity  of  the 
rebounding  shield  would  have  been  more  pronounced. 


♦ 

See  Table  2.7  for  numerical  values. 
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Tibli  2.  C 

TRIANGULAR  PUUl  (hUNLI fXAR-NOKTLASTIC ) 
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- 
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- 

4.34 

1 

4.81 
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- 
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- 
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- 
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- 

4.56 

- 
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- 

1.85 

- 

1.51 

3.01 
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1.81 

2.97 
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2.81 
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2.56 
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2.14 

- 
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THIN-BA8E  MODEL  (EXPONENTIAL  PULSE,  J4CWL1NEAK,  N0NELA6TIC) 


SECTION  3 


LINEAR  ELASTIC  ANALYSIS  — INTRODUCTION  AND  RESULTS 


3.1  Introduction 


As  mentioned  earlier,  the  purpose  of  the  analysis  presented  in 
t*'4?  and  the  following  section  is  to  provide  insight  into  the  stress- 
wave  action  which  results  in  rebound  cracking  and  spalling  and  to  pro¬ 
vide  checks  on  the  numerical  calculations  presented  in  Section  2.  As 
in  the  numerical  calculations,  the  basic  model  taken  is  a  tnin  shield 
over  a  thick  base  with  zero  bond  strength  across  the  interface.  One¬ 
dimensional  linear  wave  theory  is  used  to  calculate  the  response  of  the 
laminate  to  a  uniform  loading  with  a  step  rice  and  exponential  decay. 
Cracking  occurs  if  the  circumferential  tensile  stress  exceeds  a  critical 

value  a  ,  and  spalling  occurs  if  the  radial  tensile  stress  exceeds  a 
cr 

critical  value  0  ,  in  general  different  from  c  .  Since  the  problem 

is  linear,  the  stresses  generated  are  proportional  to  the  peek  pressure 
P.  Thus.-  for  the  circumferential  and  radial  stresses  c  fid  a  ,,  we 

V  1 

may  write 


0 


0 


?C{a  ) 


FMa  > 


where  C  and  S  depend  thiough  the  on  the  pulse  shape  factor  $ 

defined  by  p  =  Pf(t),  and  or.  material  properties.  Then  for  cracking 
we  have 


and  for  spalling 


P2P 

cr 


pap 

sp 


_®E_ 


S  v  “i  ) 


(3.1a) 


(3.1b) 
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The  relations  (3.1)  define  failure  boundaries  once  the  functions  C(a1> 
•nd  S(a^)  are  known.  The  method  for  calculating  these  functions  and 
the  corresponding  failure  boundaries  is  given  in  Section  4  for  pulses 
with  a  step  rise  and  exponential  decay.  The  remainder  of  this  section 
gives  a  comparison  of  the  failure  boundaries  so  found  with  those  from 
the  numerical  analysis  of  Section  2. 

3 . 2  Result  s 

Figure  3.1  shows  the  cracking  boundaries  for  pulses  with  a  step 
rise  and  an  exponential  decay.  The  two  upper  curves  are  from  the 
numerical  analvsls  and  are  the  same  as  those  of  Fig.  2.6.  The  lower 
curve  is  from  the  graphical  analysis  and  Is  obtained  from  Fig.  4.14, 
using  the  material  properties  and  dimensions  for  Model  l  from  Table  2.1. 


FIG.  3.1  CRACKING  BOUNDARY  FOR  EXPONENTIAL  PULSES 
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The  failure  boundary  front  the  graphical  analysis  lies  below  the 
boundary  from  the  numerical  anclysis,  with  the  difference  becoming  greater 
at  high  impulses.  The  difference  ie  attributed  to  the  effects  of  the 
load  acting  after  separation  and  to  the  contribution  of  particle  velocity 
to  the  wave  propagation  velocity,  which  art  accounted  for  in  the  numeri¬ 
cal  analysis  but  not  lr  the  graphical  analysis.  Both  of  these  effects 
would  be  expected  to  raise  the  failure  boundary  and  would  be  greater  at 
larger  impulses  (longer  pulse  duration),  In  agreement  with  Fig.  3.1. 

Figure  3.2  shows  the  graphical  spall  boundary  fot  exponential 
pulses.  The  spall  region  tc-  the  right  of  the  dashed  Hue  of  unit  slope 
is  due  to  the  effect  of  recombination  at.  the  ahield-bace  Interface.  The 
results  of  the  numerical  calculations  are  givt  .  by  the  points,  the 
crosses  indicating  spalling.  The  distribution  of  the  points,  which  was 
chosen  to  determine  the  cracking  boundary,  is  not  ideal  for  comparison 
with  the  graphical  spall  boundary.  However,  the  points  that  are  available 
agree  well.  The  most  significant  agreement  is  near  th<*  center  cf  the 
middle  square,  where  spall  and  no-spall  pointn  come  c.'ose  to  bracketing 
the  grcphic*~  Doundary.  The  difference  is  attributed  to  the  effect  of 
particle  velocity  on  wave  velocity. 

The  foregoing  comparisons  show  agreement  between  the  numerical 
and  graphical  solutions.  The  major  difference  in  the  cracking  boundaries 
from  the  two  methods  is  attributed  to  the  affects  of  the  load  acting 
after  separation  and  of  particle  velocity  on  wave  velocity.  'J  he  points 
available  from  the  numerical  analysis  a: c  consistent  with  the  graphical 
spall  boundary,  but  more  points  are  needed  to  establish  the  significance 
of  the  protrusion  attributed  to  the  effect  of  separation. 
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IMPULSE  10* - I0P» 


FIG.  3.2  SHALL  BOUNDARY  cOR  EXPONENTIAL  PULSES 


SECTION  4 


CRACKING  AND  SPALLING  BOUNDARIES  TOR  EXPONENTIAL  PULSES 
FROM  GRAPHICAL  ANALYSIS 


4. 1  Introduction 

In  tills  section  cracking  and  spalling  boundaries  are  determined 
for  exponential  pulses  using  n  aemi graphical  procedure  baaed  on  the 
method  of  characteristics. 

4. 2  One-Dimensional  Linear  Nave  Theory 

The  wave  equation  for  linear  one-dimensional  stress  waves  la 

^  .  c2i!|  <4.1, 

dt  Sx 

where  y  is  di snlacement .  t  is  time,  x  la  the  space  coordinate, 
and  c  is  the  wave  propagation  velocity.  On  the  straight  lines 
(characteristics)  in  the  x, t  plane 

x  =  ±  ct  +  constant  (4.2) 

Substituting  (4.2)  into  (4.1)  yields 

c  *  ±  4-  constant  (4.3) 

d*  dt 

Taking  p  for  pressure  and  p  for  density,  Hooke’s  law  may  be  written 

p»-  pc2  (4.4) 

ox 
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Substituting  (4,4)  snd 


u 


dv 

dt 


into  (4.3)  yields 


p  =  4  pou  +  constant  (4.5) 


Thus,  along:  the  lines  x  =  4  ct  +  constant,  the  relations  p  ®  T  peu  * 
constant  hold. 

He  no*  derive  a  basic  relation  between  loading  surface  states  at 

t  and  t  +  r,  where  t  is  t*lce  the  transit  time  of  the  shield  (2h'c). 
o  o 

Figure  4.1a  shows  an  x,t  diagram  for  a  two-t ayered  (thick  base)  laminate. 

The  characteristics  C  +  and  C  are  lines  of  x  =  4  ct  ♦  constant  in 

o  o 

+ 

the  shield,  and  the  characteristic  C.  is  a  line  of  x  =  c.  t  +  constant 

DO  D 

in  the  base.  From  the  preceding  discussion,  chsrarteri sties  in  the 

x,t  plane  map  into  characteristics  in  the  u,p  plane  (Fig.  4.1b). 

However  lr.  the  absence  of  C.  characteristics,  th**  u,p  state  behind 

b 

C.  must  be  constant.  Therefore,  the  mapping  of  C.  in  the  u.n  plane 
bo 

degenerates  to  *  point.  Moreover,  in  the  absence  of  character! sties, 

+ 

all  the  states  in  the  base  mast  lie  on  the  same  7.  •  For  *h  Initial 

D 

state  u  -  p  =  0  in  the  base,  this  must  be  p.  as  indicated.  Since 

bo 

+ 

the  states  in  the  base  are  also  interface  states,  J\  is  the  line  of 

bo 

interface  states. 

As  can  be  seen  from  Fig.  4.1a,  along  the  loading  surface  x  =  0 

the  points  t  =  t  and  t  =  t  +  t  are  connected  by  the  characteristics 
o  o 

C  and  C  ,  which  meet  at  the  interface.  Therefore,  as  shown  ic 
o  o 

Fig.  4.1b,  the  state  u  ,  p  at  t  *  t  and  the  state  u,p  at  t  =  t  + 

o  o  o  o 

i  must  be  connected  to  the  interface  state  u,p  by  T0  and  F0  • 

Hence,  we  may  write 

?  =  cbCbG 

p  “  p  *  -  pc(u  -  u  )  (4.6) 

o  o 

p  -  p  =  pc ( u  -  u) 


•JO 


o 


FlG.4.1  SCHEMATIC  *, 


u 


»  »  -  »  »«(•  -  > 


i  AND  U,  p  DIAGRAMS  WITH  NOTATION 


i 

I 

F 
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Ellmin'Atlrg  p  and  u  yields  for  the  relation  between  loading  surface 

states  at  t  and  t  *  i , 
o  o 


p  -  ecu  -  K(p  +  ycu  ) 
o  o 


(4.7) 


whei v  K  is  the  reflection  coeffic.ent  defined  as 


_  _££_ 


4.3  Exponential  Pulses 


(4.C) 


For  an  exponential  pulse, 


p  -  Pe 
o 


p  *  Po 


■X(t  +  T) 
o 


Hence, 


n  -Xt 

*—  =  e 

Po 


(4.9) 


To  simplify  the  equation-  in  the  follo^ng,  will  use 


Also,  »e  adopt  a  slope  factor  n  given  by 


(4.10) 


PQ  -  npc  uq 


(4.11) 


to  denote  the  slope  of  the  line  connecting  u0’*>o  at  t^) 

with  the  origin. 

Using  (4.9)  and  (4.11)  to  ellalnate  u  and  p  iron  (4.7)  yields 

o  o 


p  3  peu 


tie  -  (nr  1)K 


(4.12) 


For  simplicity  this  expression  Is  denoted  by  K  instead  of  Kp  as 
used  elsewhere. 
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For  given  laminate  materials,  (4.12)  gives  the  relation  between  u  and 
p  at  t  =  t  +  t,  in  terms  of  the  slope  factor  n. 

Should  separation  occur  at  t  +  0.5t,  the  intfrface  condition 
becomes  p  =  0.  This  corresponds  to  pb  =  0  and,  from  (4.8),  to 
K  =  -  1.  Thus,  if  separation  occurs,  instead  of  (4.12)  we  have 


P  = 


ne 

^°U  ne  +  (n  +  1) 


(4.13) 


The  second  equation  of  (4.6)  is  a  relationship  between  a  loading 

surface  state  u  ,p  and  the  corresponding  interface  state  u,p. 
o  o 

Using  the  first  equation  of  (4.6)  and  (4.11)  in  the  second  equation  of 
(4.6)  to  eliminate  uq  and  u  yields 


P 


K  +  1  n  +  1 
Po  2  n 


(4.14) 


For  a  given  laminate,  (4.14)  gives  the  relation  between  the  pressure  pQ 
at  the  loading  surface  at  t  and  the  pressure  p  at  the  interface 
at  tQ  +  0.5t  in  terms  of  the  slope  factor  n. 

Figure  4.2a  shows  the  x,t  diagram  for  a  two-layered  (thick  base) 

laminate  for  which  separation  does  not  occur.  In  the  region  Ola,  only 

C+  characteristics  occur.  Therefore,  each  characteristic  in  this  region 

+ 

maps  into  a  point  on  the  same  r  in  the  u,p  plane.  If  u  -  p  =  0  for 

t  <  0  (n  =  1),  this  must  be  T  +  (p  =  peu)  as  indicated  in  Fig.  4.2b. 

+ 

The  segment  of  T0  involved  is  determined  by  the  load  applied  during 
the  interval  0  £  t  <  t.  For  an  exponential  pulse,  as  shown  in  Fig.  4.2b, 
the  segment  extends  from  p  =  P  to  p  =eP.  This  gives  the  loading 
surface  states  for  0  £  t  <  t. 

To  determine  the  interface  states  for  0.5t  £  t  £  1.5t,  we  use 
(4.14)  with  i\  =  1  and  obtain 

p(t)  =  (K  +  l)Pe(t  -  0.5T)  (4.15) 

These  states  lie  on  the  line  p  =  Pbcbu  along  the  segment  12  ,  as  shown 
in  Fig.  4.2b. 
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•e-un-i 

FIG.  4.2  x,t  AND  u,p  DIAGRAMS  FOR  EXPONENTIAL  PULSES 
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The  loading  surface  states  for 


t  4  t  <  2t  are  found  using  (4.12) 

with  n  a  J  (since  u  ,p  lies  01,  p  =  pcu).  This  gives 

o  o 


p  = 


DC  U 


(4.16) 


Thus,  these  states  lie  on  a  straight  line  of  slope  oc[e/(e  -  2K)j  passing 
through  the  origin. 


The  lines  I  and  II  o'.  Fig.  4.2b  show  two  possible  locations  of 
this  line  with  i-espect  to  the  line  p  =■  -  pcu.  To  facilitate  the 
discussion  oi  these  case;,,  xe  first  obtain  the  equations  for  the  loci  of 
subsequent  loading  surfp.ee  states. 


To  obtain  the  line  of  loading  surface  states  for  27  4  t  <  3t  ,  we 
put  n  -  e/(e  -  2K )  =  \  froti  (4.16)  into  (4.12)  and  obtain 


P  ~ 


pcu 


_ xe 

xe  -  (x  +  1)K 


(4.17) 


by  continuing  this  process  we  obtain  for  the  successive  loading  surface 
states 


where 


p  = 

pcu 

P  2  p  eP 

p  = 

xpeu 

2 

eP  2  p  >  e  P 

p  = 

yipcu 

2  3 

e  P  2  p  >  e  P 

p  - 

y2rcu 

e^P  2  p  >  e^P 

etc , 

1 

X  3 

1-2* 

X 

y\  x  -  (1  +  xK 

_  ~'l _ 

y2  "  yx  -  (1  +  yjl? 


etc . 


14. ie) 


(4. 19) 


and 


•=  K/e. 
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As  shown  bel«-w,  if  separation  occurs  it  takes  place  first  at 
t  i  1.5t.  Then,  U8i;*g  14.13)  instead  of  (4.12),  we  obtain  tor  the  se- 
cuence  of  loading  surface  states 


»hc-re 


p  =  pcu 

P  =  xocu 

P  =  z^ccu 

p  =  ZjC-cti 

etc . 


P  2  p  >  eF 

2 

eP  2  p  >  e  P 

2  3 
e  P  2  p  >  e  P 

3  .  4 

e  P  £  p  >  e  P 


_ xe _ 

Z1  xe  +  <1  +  x) 


*2  V  -  (1  +  *j> 

etc . 


(4. 2C) 


(4.21) 


4 •  4  txponentiai  ^itsea  with  Slow  L'atv  (h  u) 

It  Is  readily  shown  that  for  X  <  e  (corresponding  to  I  in 
Fig.  4.2b>  tensile  stresses  do  not  develop  in  the  shield  or  base  *nd 
hence  cracking  or  spalling  cannot  occur.  The  u,p  diagram  for  K  <  e, 
extending  beyond  that  of  Fig.  4.2b.  is  shown  in  Fig.  4.3;  the  x,t 
diagram  is  the  same  as  that  of  Fig.  4.2a,  Consac.  tive  loci  of  u,p  states 
at  the  loading  surface  are  the  line  segments  on,  ub ,  be,  etc.  These 
correspond  to  the  states  given  by  (4.18). 

The  seijucuCc  of  lines  ot,  cb ,  c“ ,  etc .  cf  Fig.  4.3  "ppears  to 
approach  a  limit.  The  slope  of  this  line  is  easily  found  from 
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u 

oc  m* -  * 

FIG.  4.3  v,p  DIAGRAM  FOR  A  PULSE  WITH  LOW  RATE  OF  DECAY 

(4.19)  by  putting  yn+J  =  yn-  This  yields  for  the  equation  of  the  limit¬ 
ing  line, 


p  =  U-S  ecu  (4.22) 

e  —  K 

Hence,  as  can  be  seen  from  Fig.  4.3,  for  K  <  e  all  the  u,p  states  in 
the  shield  and  the  base  lie  within  the  region  bounded  by  the  lines 
p  ~  qcu  and  p  —  —  pcu.  Hence  no  tension  develops  in  the  shield  cr 
base. 

4.5  Exponential  Pulses  with  K  =  e 

The  x,t  and  u,p  diagrams  for  the  special  case  K  =  c  are  given 
in  Fig-  4.4.  The  triangles  of  state  points  a2b,  c**d,  etc.,  of  Fig.  4.3 
have  all  become  segments  of  the  line  p  =  -  pen  and  the  triangles  of 
state  points  b3c ,  d5e,  etc.,  have  become  segments  of  the  line  p  =  pcu. 

All  state  points  for  the  triangular  areas  2b3 . are  represented  by 

the  origin  u  =  p  =  0.  Since  pressures  are  always  positive,  cracking 
or  spalling  does  not  occur.  For  K  -  e,  the  limiting  line  (4.22) 
approaches  the  p-axis. 
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AND  u,f  DIAGRAMS  FOR  K 
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4.6  Exponent!*!  Pulsis  with  K  >  e 


As  shown  above,  for  K  <  e  no  tensile  stresses  develop  In  the 
shield  or  base.  For  K  >  e,  tensile  stresses  always  occur,  as  can  be 
readily  shown  trom  Fig.  4.2b. 

In  Fig.  4.2b,  the  locus  of  interface  Staten  for  1.5t  i  t  <  2.5t 

in  the  segment  of  p  =  C.c,u  (line  of  interlace  states)  which  would  be 

b  b 

intersected  by  the  extension  of  the  lines  emanating  from  II  Hence, 
the  first  possibility  of  tensile  stress  occurs  at  t  =  1.5t.  For  a 
perfect  bond,  the  magnitude  cf  the  tensile  stress  at  t  -  1.5t  is 
given  by  the  intersection  of  the  line  through  a2  with  the  line  of  inter¬ 
face  states.  For  the  interval  1.5t  i  t  <  2.5t,  the  tensile  stress  with 
a  perfect  bond  is  t  maximum  for  t  =  1.5r. 

For  ze-r  bond  strength,  tensile  stresses  at  the  interface  cannot 
be  sustained  and  separation  occurs  at.  t  =  1.5t.  Subsequent  interface 
states  lie  on  the  segment  of  p  =  0  intersected  by  the  lines  emanating 
f  rom  1 1 , 

For  zero  bond  strength,  the  region  K  >  e  divides  Into  two  parts 
for  which  the  u,p  diagrams  are  basically  different.  If  point  2  of 
Fig.  4.2u  lies  above  p  =  0,  as  shown,  one  type  of  diagram  results;  if 
it.  lies  below  p  -  0,  a  different  type  results.  The  division  occurs  for 
point  2  lying  on  p  =  0.  Puic‘lc  for  which  point  2  lies  below  o  -  0 
are  pulses  with  fast  decay  ;  pulses  for  which  point  2  lies  above  p  =  0 
are  pulses  with  moderate  decay. 

The  x,t  and  u,p  diagrams  for  p2  =  0  are  shown  in  Fig.  4.5.  As 
can  be  seen  ir  the  figure,  state  2  is  determined  by  the  intersection  of 
the  lines 

p2  -  p?  =  cc(u2  •  y 

P2  ~  P_  =  "  P=<u2  "  u_  '> 

a  a 
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Xp.CU 


er 


P_  = 

a 


a 


with  x  from  (4.21),  and  p^  =  0  we  obtain  the  condition 


K 


(4.23) 


As  Indicated  in  Fig.  4.5b,  the  tinal  loading  surface  state  is 
U  -  p  =  0.  Hence  the  shield  separates  and  is  brought  to  rest  by  the 
pulse.  As  is  shown  below,  for  pulses  with  high  decay  rate  the  shield 
is  left  with  a  final  velocity,  and  for  pulses  vi+h  moderate  decay,  the 
shield  is  driver,  back  toward  the  hast-  and  recombination  occurs.  The 
regions  in  the  K,e  plane  for  slow,  moderate,  and  fast  decay  are  shown 
later  in  Fig.  4.15. 

4, 7  Exponential  Pulses  with  Fast  Decay  (0  <  e  <  K/K  *  1 ) 

Typical  u,p  diagrams  for  pulses  with  fast  decay  are  shown  in 
Fig.  4.6;  the  x,t  diagram  is  the  same  as  in  Fig.  4.5a.  Figure  4.ba  gives 
the  overall  aspects  of  the  diagram  and  Figs.  4.6b  and  c  show  details  for 
relatively  large  and  small  values  of  e.  It  is  readily  seen  from 
Figs.  4.6b  and  c  that  the  sequence  of  triangular  regions  a2b,  b3c,  etc. 
reaches  the  limiting  region  V.4N,  and  the  sequence  2'o3,  3c4,  etc.  reaches 
Nn.\ .  These  two  limiting  regions  give  the  final  state  of  notion  and 
stress  in  the  shield.  As  can  be  seen  from  the  x,t  diagram  of  Fig.  4.5a, 
.fter  t  =  1.5t  the  motion  is  periodic  with  period  t.  Since  velocities 
are  negative  throughout,  the  shield  remains  separated. 


The  maximum  tensile  stress  (a  =  P„)  indicated  in  Fig.  4.6  is 

max  z 

found  from  the  equations 


p2  “  Pg  "  Pp(p2'  “a1 

p 2  -  Pa  =  PC<V 

-  a 

p2  *  Pbchu2  *  (1  *  K)eP 


p_  =  xpcu_  »  eP 
a  a 


and  is  given  by 


-P2  ■  ■  [K  -  <1  hk  K>«]p 


(4.24) 


(4.25) 


For  the  special  case  (1  +  K)e  =  K  of  Fig.  4.5,  this  reduces  to  a  "  0 

max 

as  it  should. 


For  the  limiting  case  e  -♦  0  (very  short  pulse)  we  find  a  =  KP. 

max 

The  corresponding  u,p  diagram  is  shown  in  Fig.  4.7.  The  regions  a2b,  b3c, 
etc.,  all  coincide  with  the  region  NNN,  and  2b3,  3c4,  etc.  coincide  with 
NnN.  The  diagram  of  Fig.  4.7  applies  to  short  pulses  of  arbitrary  shape, 
i.e.  all  pulses  with  step  rise,  a  continuous  decay,  and  a  total  duration 
T  such  that  T  £  t. 

4.8  Exponential  Pulses  with  Moderate  Decay  [k/(K  +  1)  <  e  ■<  k3 

The  analysis  of  pul Be s  with  moderate  decay  is  more  complicated 
than  that  for  pulses  with  slow  or  fast  decay  because  recombination 
occurs.  Figure  4.8  shows  the  u,p  diagram  for  a  pulse  with  e  slightly 
less  than  K.  The  x,t  diagram  (until  recombination  occurs)  is  the  same 
as  that  of  Fig.  4.5a. 
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FIG.  4.7  U;P  DIAGRAM  FOR  e  -  0 


The  regions  of  states  ola,  la2,  and  a2b  are  found  as  before. 

Since  e  <  K,  separation  of  shield  and  base  occurs  at  t  =  1.5t.  Hence 
the  region  of  states  2b3  is  determined  by  the  condition  p  =  0  at 
interface.  The  following  two  regions  b3c  and  3c4  are  also  easily  de¬ 
termined. 

From  the  x,t  diagram  of  Fig.  4.5a  it  is  seen  that  the  segment  of 
the  u-axis  23  in  Fig.  4.8  gives  the  interface  velocity  for  1.5t  £  t  < 
2.5t,  and  the  segment  34  gives  the  interface  velocity  for  2.5t  ^  t  <  3 
Since  the  velocities  corresponding  to  the  segment  23  are  smaller  than 
those  corresponding  to  the  segment  34,  the  outward  displacement  of  the 
shield  during  1.5t  £  t  s  2.5t  is  less  than  the  possible  inward  dis¬ 
placement  during  2.5t  st  £  3.5t.  Hence  recombination  of  the  base 
occurs  at  an  instant  t  =  cpT,  where  2.5  <  cp  <  3.5. 
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FIG.  4.8  u,p  DIAGRAM  FOR  A  PULSE  WITH  MODERATE  DECAY  K  '  «  >  e , 
e  [2  •  e)  ■*  2K.  Recombination  ol  t  <  3.5r 

Alter  the  instant  t  =  rjr ,  state  points  for  elements  of  the  shield 
at  the  interface  are  no  longer  determined  bv  the  condition  p  ~  0:  in¬ 
stead  they  must  lie,  as  before  separation,  on  the  line  p  =  s^c^u.  Hence, 
the  triangle  ol  state  points  3c4  is  not  completely  realized;  only  points 
within  the  legion  33 'c 'c  actually  occur.  Tie  remaining  part  3 'c '4  of  the 
region  3c4  is  transferred  to  the  interface  line  p  =  a8  Indicated 

in  Fig.  4.8. 

A  more  detailed  version  of  the  diagram  ol  Fig.  4.8  including  later 
states  is  shown  In  Fig.  4.9.  Inspection  of  the  interface  states  indi¬ 
cates  that  separation  and  recombination  occur  repeatedly. 
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Figures  4,8  and  4,9  are  valid  for  initial  recombination  during  the 
interval  2,5t  <  t  <  3.5t.  Figure  4,10  gives  the  u,p  diagram  for  re¬ 
combination  at  t  =  3.5t,  This  occurs  when  the  line  of  states  be  coin¬ 
cides  with  the  p-axis,  since  then  the  velocities  corresponding  to  the 
line  segments  23  and  34  on  the  u-axis  are  equal  and  opposite.  For  the 
line  of  states  be  to  coincide  with  the  p-axis,  the  slope  factor  z of 
(4.20)  must  be  infinite.  This  gives 

2K  =  e(e  +  2)  (4.26) 

The  parameters  K  and  e  of  Fig,  4.10,  in  addition  to  satisfying 
the  condition  z  -♦  *  ,  were  chosen  so  that  after  recombination,  i,e,, 
subsequent  to  the  region  c4d,  the  u,p  diagram  is  identical  to  that  of 
Fig.  4.5b.  This  occurs  wnen  state  5  (corresponding  to  state  2  of 
Fig.  4.5b)  is  on  the  u-axis.  This  results  in  the  condition 

e3  +  |e2  +  |  =  1  (4.27) 


FIG.  4.10  u,p  DIAGRAM  FOR  THE  SPECIAL  CASE  £[l  +  (e/2)]  =  K, 
e3  +  (3/2)*2  +  e/2  -  1  (e  -  0.6,  K  -  0.5  S) 
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This  equation  has  one  real  root  in  the  interval  0  <  e  <  1  which  is 
approximately  0,6,  The  corresponding  value  of  K  frojn  (4.26)  is  0,78. 
Thus,  for  e  =  0,0,  K  =  0,78,  the  diagrams  of  Fig.  4,5b  and  that  of 
Fig.  4.10  become  identical  subsequent  to  region  c4d.  This  means  that 
the  shield  eventually  comes  to  rest. 

Figure  4.11  gives  the  u,p  diagram  for  which  initial  recombination 
occurs  during  the  interval  3.5t  <  t  <  4.5t.  For  this  case  recombination 
does  not  occur  prior  to  t  =  3.5t,  because  the  average  outward  velocity 
corresponding  to  the  segment  23  on  the  u-axis  is  greater  than  the  average 
inward  velocity  corresponding  to  the  segment  34.  However,  the  average 
velocity  of  the  two  segments  34,  45  is  greater  than  half  the  average 
velocity  of  23,  hence  recombination  occurs  for  t  <  4.5t.  The  regions 
4d5  and  d5s,  respectively,  in  Fig.  4.11  correspond  to  the  regions  3c4 
and  c4d  of  Fig.  4.8  and  are  not  fully  realized. 


FIG.  4.11  u,p  DIAGRAM  FOR  A  PULSE  WITH  MODERATE  DECAY. 
Recombination  at  t  >  3.5r. 
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I:.  Fig.  4.9  tho  magnitude  of 
In  the  ah. eld  for  the  states  shown 
ponding  to  the  states  y  and  z. 
subsequent  tensile  stresses  should 
and  /.  are  taken  as  bounds  on  the 


the  maximum  tensile  stl^ss  generated 
lies  between  the  stresses  corres- 
81  nee  the  load  decreases  continually, 
be  less.  Hence,  the  stresses  at  y 
maximum  tensile  stress  generated. 


State?;  y  and  z  are  determined  by  the  intersection  of  the  line 
xs  in  Fig.  4.8  wihh  the  lines  p  =  ecu  and  p  =  •  The  equation 

of  the  line  xs  te 


P  -  Px  =  -  pc(u  -  ux> 


with 


=  e4P  =  mpeu^ 

and  m  given  by  (4,12)  foi  n  -  from  (4,21).  These  equations 
yield 

-  py  =  cy  "  Pl(e3  4  e2  +  e  -  K)  K  -•  e4j  (4.28) 


and 


y..  - 


The  stresses  given  by  (4,28)  and  (4.29)  are  valid  for  initial  recombina¬ 
tion  occurring  in  the  interval  2.5t  it  <  3, St. 


From  Fig.  4.11,  for  initial  recombination  In  the  interval  3. St  * 
t  <  4.5t  ,  the  stresses  at  y  and  z  srr  found  from  the  equations 


P  -  Px  “  -  pc(u  -  ux> 


wl  th 


*  e  P  »  mfeu 


SI 


Iron  (4.21  >. 


These  yield 


F 

I 


\, 

m* 

l! 


r 


and  m  given  by  (4.12)  foi  n  = 


-  p  -  0  =  pice*  ie^  +  e-  K)K-  e*] 

y  y 

and  the  expression  of  (4.29)  for  a  . 

Equations  (4.28)  and  (4.29)  suggest  the  general  form 

iTi  r i  —  Z  t'  -  2  3  2 

=  PI  (e  *  e  +  e  •••  .  .  .  .  ♦  e  +  e  +  e  -  K  )K 

y 


(4.30) 


*] 


(4.31) 


.  here 


•  N  +  1 


for  Initial  recombination  occurring  in  the  Interval 

(N  +  0.5)-  St  S  (N  +  1. 5)t  V.  =  2,3,4 . 


For  large  N  ,  a  slmpl.  r  ar.fl  more  accurate  estimate  of  the  maxi¬ 
mum  tensile  stress  can  be  obtained  as  followt .  The  state  points  a,  b, 
c,  etc.  of  Figs.  4.8  and  4.11  are  connected  by  a  straight  line  ol 
slope  -(1  -  e)/(l  +  e),  as  are  the  points  b,  c,  d,  etc.  For  the  first 

set  cf  these  points  this  can  oe  readily  veriiied  b_»  uull» K  that  they 

2  2  3  3 

have  the  coordinates  (ep/xpc,  eP),  (e  P/7  pc,  e  P),  (e  P/7  pc,  e  P), 

4  4  *  * 

(e  P/7,jpc,  e  P),  etc.  In  the  u,p  plane,  with  x,  z^,  z2>  Zy  etc.  as 

given  by  (4.21).  Hence,  the  equation  for  the  straight  line  ab  is 


P 


1  - 
1  r 


peu 


e  -  (1  -  e)K 
1  +  e 


and  is  satisfied  by  the  coordinates  of  the  points,  c,  d,  etc  .  In 
Fig.  4.11,  intersection  ol  this  lina  with  the  u-axls  yields  the  point 
(ii,0>,  vi  th 


peu 


o  - 


(1  ~  e)K 
1  -  e 


2P 


(ItK)e-K  2p 
1  -  e 
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For  decreasing  values  of  e  the  intersection  point  (u,0)  moves  toward 
the  origin,  reaching  the  origin  for  e  =  K/(l  +  K).  As  u  approaches 
the  origin,  more  velocity  increments  represented  by  line  segments  of  the 
«et  34,  45,  etc.  are  needed  to  nullify  the  outward  velocity  corresponding 
to  the  line  segment  23,  implying  that  recombination  will  be  continuously 
delayed.  Moreover,  ar  u  approaches  the  origin,  the  region  Uwxs  (Fig. 4. 11) 
approaches  the  limiting  region  uw'x's'.  Hence,  for  e-*K/(l  +  K)  the 
maximum  stress  generated  in  the  shield  is  bracketed  between  the  values 
C  ,  and  o  t 

y  * 

For  a  <  we  have 
z 


0  ,  =  (1  +  K)  c  i 

z  y 


The  intersection  of  the  line  x'a'  with  p  =  pc(u  -  u)  through  u,0 
gives  p  t  =  -  a  /  .  The  line  x's'  is  found  using 

y  y 


p  =  -  pc ( u  -  ux) 


r>  ,  =  pc(u,-u,)  =  -pc(u/-u) 

W  W  X  w 


P  r  =  c.  C.  U  , 
w  b  b  w 


which  yield 


P 


pcu 


(1  +  K)e  -  K 
1  -  e 


2KP 


The  intersection  of  this  line  with  p  =  pcu  gives 


"  P. 


J  f 

y 


(1  -*  K )e  -  K 
1  -  e 


KP 


(4.32) 


For  e  ■  K/(K  +  1),  (4.33)  gives  ,  =  0.  This  is  as  expected  from 
the  u,p  diagram  for  this  case  in  Fig.  4.5. 


i 


i 
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Equations  (4. 31).  (4.32),  and  (4.29)  provide  estimates  ol  the 
maximum  radial  tensile  stress  for  pulses  ol  moderate  decay.  In  Section 
4.10  these  equations  are  used  to  calculate  the  spall  bounaary. 

4. 9  Cracking  Boundary  for  Exponential  Pulses 

To  determine  the  cracking  boundary,  we  regard  the  shield  as  a  ring 
in  uniform  radial  motion.  Cracking  occurs  if  the  outward  displacement 
is  sufficient  to  develop  tensile  stresses  exceeding  the  fracture  stress 
of  the  shield.  For  pulses  with  slow  decay  (K  s  e),  .separation  does  not 
occur;  hence,  cracking  Is  not  possible.  However,  for  pulses  with  moder¬ 
ate  or  last  decay  (  K  >  e),  separation  always  occurs  and  cracking  is 
possible  if  the  load  amplitude  is  sufficient. 

To  determine  the  maximum  outward  displacement,  tic  average  out¬ 
ward  velocity  of  the  shield  at  the  separation  preceding  cracking  is 
required.  For  pulses  with  fast  decay,  permanent  separation  occurs  at 
t  =  1.5t.  For  pulses  with  moderate  decay,  the  average  outward  velocity 
of  the  shield  is  greatest  at  first  separation,  i,e.,  at  t  =  1.5t.  This 
car.  be  seen  in  the  example  of  Fig.  4.9  by  comparison  of  the  average 
velocities  at  and  (the  velocities  av.  s^  fall  within  the  left¬ 

most  triangle)  Since  separation  velocity  is  proportional  to  the  load 
amplitude,  and  the  effective  amplitude  decreases  for  each  successive 
separation,  the  average  outward  velocity  will  be  greatest  at  first 
separation,  i.e.,  at  t  =  1.5t.  Therefore,  the  average  velocity  needed 
to  determine  the  maximum  outward  displacement  is  that  at  t  =  l.Sr  foj 
all  pulses. 

The  average  velocity  at  t  =  1.5t  is  readily  found.  The  velocity 
of  the  shield  at  the  interface  at  t  =  1.5'  is  th  it  at  state  2  in 
Fig.  4.8.  The  velocity  at  the  loading  surface  av  the  same  Instant  is 
on  the  line  ah  of  Fig.  4.8  and  is  given  by 

P(1.5t)  =  =  xpeu^  ■  (4.33) 

with  x  from  (4.19).  Since  for  exponential  pulses  the  transformation 
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x,t  -*  u,p  is  linear,  for  points  between  the  loading  surface  and  the 
interface  the  velocity  will  lie  on  the  straight  line  in  the  u,p  plane 
connecting  the  loading  surface  state  and  the  interface  state.  Thus, 
the  average  velocity  of  the  shield  is  simply 

v  =  \  <u2  +  V>  (4.34) 


From  Fig.  4.3  we  have 


Together  with  x 


P_  '  P2  =  "  Oc(u_  -  u2) 

a  a 

P?  "  p2  =  CC(U2  ’  U2) 

p2  =  (1  +  K)p^  =  pbcbu2 

p_  =  pa  =  eP 
a 

=  xqcu_ 

a  a 

from  (4.19)  these  equations  >ield 
pcu2  =  t(l  -  K)fc  -  XjP 


(4.35) 


(4.36) 


- V.  <  -i  ^  —  t  «  Oi?\  +  t  /  A  oc  \ 

V/OdlMa  V  ■  •  VV  >  n  *  VM  \  •  »  >'*'  /  6*  •  W.U 


-  A  —  [k  -  (1  -  K)e  +  (2K  -  e)e1/2J  (4.37) 

2  pc 

The  outward  displacement  of  the  shield  is  determined  by  the  average 
outward  velocity  at  separation  and  the  subsequent  load.  Here  we  neglect 
the  subsequent  load  and  thus  obtain  a  lower  bound  on  the  load  required 
for  cracking.*  The  kinetic  energy  (per  unit  area)  of  the  shield  st 


To  account  for  the  subsequent  load  is  tif  difficult,  out  it  »a»  not 
accomplished  in  tine  for  inclusion  in  this  report. 
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separation  is 


Kinetic  energy  =  phv2 


where  h  is  the  shield  thickness.  For  a  ring  in  uniform  radial  motion, 


the  strain  energy  is 


1  h  2 

Strain  energy  =  -  —  n 

“  E9  9 

where  a.  is  the  hoop  stress  and  E  is  the  corresponding  elastic 

0  9 

modulus.  Equating  the  kinetic  energy  to  the  strain  energy  gives 


2  -  Zi 

V  pE 


2  2 

9  p  c9 


(4.38) 


where  c  is  sound  speed.  For  a  homogeneous  isotropic  material  the 

9 

ratio  of  the  radial  (dilatational)  sound  speed  c  to  the  circumferential 

(plate)  sound  speed  ca  is 

d 


g  _  _£_  =  _ 3  ~  v 

C0  VI  -  2v 


(4.39) 


Introducing  (4.39)  into  (4.38)  yields 


v  =  0Q 
pc  0 


(4.40) 


Using  (4.40)  in  (4.37),  and  putting  o.  =  O  ,  gives 

0  cr 


=  23[k  -  <1  -  K)e  +  (2K  - 


a 


Oe1/2]'1 


cr 


(4.41) 


Comparing  (4.41)  with  (3.1a),  we  have 


C  =  ~[K  -  (1  -  K)e  +  (2K  -  e)e1/2] 


(4.42) 


v  of  (4.40)  is  equal  to  -  v  of  (4.37). 
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The  cracking  boundary  (4.41)  is  shown  in  Fig.  4.12  for  K  =  0.7  and 
0  =  1 . 3 ( v  =  0.4). 


FIG.  4.12  CRACKING  BOUNDARY  FOR  EXPONENTIAL 
PULSES  IN  TERMS  OF  e 


For  an  exponential  pulse  the  impulse  1  is  given  by 

I  =  Z  =  py  •  (4.43) 

A, 

where  T  As  the  characteristic  time  duration  (see  Fig.  2.3).  Using 
■In  e  =  -  \t,  we  have 


T  _  2 

t/2  Ln(l/e) 


(4.44) 


The  cracking  boundary  in  terms  of  T/(t/2)  is  given  in  Fie.  4.13. 
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FIG.  4.13  CRACKING  BOUNDARY  FOR  EXPONENTIAL  PULSES 
IN  TERMS  OF  T 


To  obtain  the  cracking  boundary  in  the  P,I  plane  we  use  the 
transformation 


P  T 

T  -  p  T  =  cr 

cr  cr  ^n(l/e) 


(4.45) 


obtained  from  (4,43)  and  (4.44).  The  limiting  value  of  I  ,,  for 

cr 

e  -*  0  (very  short  pulses)  is  raauily  found.  For  pulses  of  arbitrary 
shape  with  durations  less  than  t,  the  momentum  at  separation  is  the 
reflected  Impulse  I  given  by 


phv 
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Dropping  the  m: nus  sign  and  putting  h  = 


c  -  gives 


cr 


1  T 

=  K  pCV  2 


Substituting  v  from  (4.38)  yields 


cr 


1  c  T 

K  c  Ccr  2 

0 


With  0  =  c/ca  and  I  =  g/K  we  obtain 

U 


I  =fa  -  =  i*a  -~ 

cr  K  cr  2  cr  cr  2 


(4.46) 


for  the  limiting  impulse  of  the  cracking  boundary  as  e  -»  0. 


Using  (4.45)  and  (4.46)  yields  the  cracking  boundary  in  the  P,I 
plane  shown  in  Fig.  4.14,  where 

/ 

*  I 

1  =  — -  (4.47) 

Vi 

The  boundary  has  the  expected  shape.  In  Fig.  3.1,  the  cracking  boundary 
of  Fig.  4.14  is  compared  to  the  cracking  boundary  obtained  numerically. 

4. 10  Spall  Boundary  for  Exponential  Pulses 

Figure  4.15  shows  the  regions  of  the  K,e  plane  for  exponential 

pulses  with  slow,  moderate,  and  fast  decay.  For  pulses  with  slow 

decay  (K  £  e)  radial  tensile  stresses  do  not  occur;  hence  spall  is  not 

possible.  For  pulses  with  moderate  or  fast  decay  (K  >  e)  radial  tensile 

stresses  always  occur  and  spall  is  possible  if  the  load  amplitude  is 

sufficient.  To  facilitate  the  discussion,  we  regard  K  as  fixed,  say 

K  =  0.7.  Then,  to  establish  the  spall  boundary,  P  /a  must  be 

*  sp  sp 

determined  along  the  line  K  =  0.7  in  Fig.  4.15. 
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FIG.  4.14  CRACKING  BOUND, 
IN  PRESSURE-IMP1 


10 


FIG.  4.15  SLOW,  MODERATE,  AND  FAST  DECAY  REGIONS 
IN  K,e  PLANE 


The  segment  of  the  line  K  =  0. 7  extending  upward  from  e  =  K 
corresponds  to  pulses  of  slow  decay,  for  which  spall  does  not  occur. 
Using  the  transformation  of  (4.44)  with  e  *  K  gives  the  vertical 
boundary  of  Fig.  4.16.  Since  larger  values  of  e  correspond  to  larger 
values  of  T,  spall  cannot  occur  in  the  region  to  the  right  of  the 
vertical  boundary. 

The  segment  of  the  line  K  =  0. 7  extending  downward  from 

e  =  K(K  +  1)  in  Fig.  4.15  corresponds  to  pulses  with  fast  decay.  For 

this  segment  we  have  from  (4.25),  with  o  -  a  .  P  =  P  , 

max  sp  sp 

P  , 

■~E  =  [g  -  (K  +  De}”  (4.48) 

»P 

Hence,  from  (3.1b), 

S  =  K  -  (K  +  l)e  (4.49) 
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FIG.  4.16  SPALLING  BOUNDARY  FOR  EXPONENTIAL  PULSES 
IN  TERMS  OF  T 


Equation  (4.49)  together  with  the  transformation  of  (4.44)  yields  the 

fast-decay  boundary  of^Fig.  4.16.  For  the  limiting  value  e  -»  K/(K  +  1), 

(4.49)  yields  P  -»  «,  corresponding  to  the  verti^  1  asymptote  in 
sp 

Fig.  4.16.  For  e  -*  0,  (4.49)  gives  P  /c  =  1/X  =  1.47,  corresponding 

sp  sp 

to  the  intersection  of  the  spall  boundary  with  the  vertical  axis  in 
Fig.  4.16. 

The  segment  of  the  line  K  =  0.7  between  e  =  K  and  e  =  K/(K  +  1) 
corresponds  to  pulses  of  moderate  decay.  For  the  part  of  this  segment 
X  >  e  >  e  ,  with  e(e  +  2)  =  2K  from  (4.26),  the  spall  boundary  is 
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and  N  =  2, 


/ 

found  using  (4.29)  and  (4.31)  with  a  =  O'  ,  P  =  P  , 

z  sp  sp 

corresponding  to  recombination  in  the  interval  1.5t  £  t  <  2.5t.  Making 
these  substitutions,  we  have 

£  (e3  +  e2  +  e  -  K)K  -  e4]"1  (4.50) 

U  K  '  £ 

sp 

Together  with  the  transformation  (4.44),  (4.50)  gives  the  spall  boundary 
for  K  >  e  >  e. 


To  determine  the  intersection  of  the  moderate-decay  boundary  with 
the  slow-decay  boundary  in  Fig.  4.16,  we  put  e  =  K(1  -  6),  6  <  <  1, 
into  (4.50)  and  obtain 


lim  =  - - - 5  (4.51) 

6  -  0  flP  (K  +  1)K 

For  K  =  0.7  this  gives  P  /a  =1.71  at  the  intersection  of  the 

sp  sp 

moderate-decay  and  slow-decay  boundaries. 


To  determine  P  /a  from  (4.50)  at  e 
sp  sp 

e  <  <  1,  and  obtain 


e  ,  we  put  e  =  e(l  +  e), 


lim 
e  -+  0 


1 

K  +  1 


+  e 


r2vi-l 


(4.52) 


For  K  =  0.7,  e  =  0.55,  and  (4.52)  yields  P  /a  =  4.5.  Thus,  as  can 

sp  sp  _ 

be  seen  from  Fig.  4.16  the  spall  boundary  for  K>e>e,  i.e.,  N  =  1, 
extends  from  the  slow-decay  boundary  to  T/(t/2)  =  3. 


✓ 

For  the  next  Interval,  N  =  3,  (4,29)  and  (4.31)  yield,  with  a  = 

z 

and  P  =  p  , 
sp 


Psp  1 

a  =  R  +  i 
sp 


r.  4  ^  3  2 

I  (e  +  e  +  e  + 


-  K)K  -  e  ] 


5i-l 


(4.53) 
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For  e  -»  e,  (4.53)  gives 


Psj^  1  re3  .  5  -  -2  -3 i-l 

a  =  iTTx  (1  +  5  e  +  c  +  c  ] 

sp 


(4.54) 


For  subsequent  intervals,  we  use  (4.29)  and  (4.32),  with  a  =  c 

z  sp 

and  P  =  P^,  to  determine  the  spall  boundary.  This  gives 


ap  _  1  (K  +  l)e  -  K 

O  K+l  1  -  e 

sp 


(4.55) 


The  spall  boundary  in  the  P,I  plane  is  obtained  from  Fig.  4.16 
using  the  transformation  (4.45)  written  as 


P  T 

i  =  P  l  a  -  .seL, 

sp  sp  2  ind/e) 


(4.56) 


together  with  the  definition 


* 

I  = 


l7cr  i 


(4.57) 


The  spall  boundary  in  the  P,I  plane  is  shown  in  Fig.  4.17.  The  spall 
boundary  of  Fig.  4.17  is  compared  with  numerical  calculations  in  Fig.  3.2, 
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SECTION  5 


EXPERIMENTAL  VERIFICATION  OF  REBOUND  CRACKING  THEORY 

t . 1  Introduction 

To  attempt  to  /erily  the  theory  of  rebound  cracking,  an  experi- 

♦ 

mental  program  was  undertaken  using  a  shield  of  phenolic  setting  resin 
and  a  thick  base  of  steel.  The  main  element  of  the  experiment  is  shown 
in  Fig.  5.1.  On  the  central  steel  mandrel  are  mounted  a  phenolic  test 
section,  two  segmented  aluminum  rings,  a  segmented  phenolic  ring,  and 
two  end  spacers.  The  segments  are  shown  more  clearly  in  Fig.  5.2. 


(a)  OUTFITTED  MANDREL  (b)  HE  PLACEMENT 


FiG.  5.1  MAIN  ELEMENT  OF  EXPERIMENT  TO  VERIFY  REBOUND  CRACKING  THEORY 


Micarta,  Westinghouse. 
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FIG.  5.2  SEGMENTED  RINGS 

To  keep  the  outer  surface  unobstructed  for  photographic  observations, 
the  explosive  load  was  applied  to  the  inner  surface  of  the  hollow  steel 
mandrel.  The  main  purpose  of  the  experiment  was  to  determine  (1)  the 
ability  o l  the  theory  to  predict  the  velocity  of  the  phenolic  segments 
and  (2)  the  conditions  under  which  the  test  section  would  crack. 

Selected  frames  from  a  framing  camera  sequence  of  a  typical  shot 
are  shown  in  Fig.  5.3.  The  segment  velocities  from  such  photographs 
were  used  to  check  the  predictability  of  the  velocity  of  the  phenolic 
segments  and  to  determine  the  input  pulse  to  the  test  section.  To 
check  the  predictability  of  the  velocity  of  the  phenolic  segments,  the 
wave  transit  time  of  the  thinner  aluminum  ring  was  made  equal  to  that 
of  the  phenolic  ring;  as  shown  below,  this  results  in  the  same  input 
pulse  to  both  rings.  Using  one-dimensional  linear  wave  theory,  the 
phenolic  velocity  predicted  from  the  velocity  of  the  aluminum  segment 
is  found  to  agree  with  the  observed  segment  velocity  to  within  5  percent. 
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The  initial  velocity  of  the  test  section  is  taken  to  be  the  same 
as  the  velocity  of  the  phenolic  segment,  which  is  of  equal  thickness. 

It  is  found  that  the  predicted  cracking  velocity,  based  on  the  static 
stress-strain  curve  for  the  phenolic,  and  observed  cracking  velocity 
agree  to  within  10  percent.  A  cracked  test  section  is  shown  in  Fig.  5. 4. 

The  velocity  of  the  thicker  aluminum  ring,  together  with  that  of 
the  thinner  ring,  is  used  to  obtain  information  of  the  shape  of  the 
input  pul se  at  the  segment-mandrel  interface.  It  is  found  that,  of  the 
triangular  and  exponential  shapes,  the  pulse  is  matched  best  by  a  trian¬ 
gular  shape. 

Additional  information  on  the  experimental  setup  and  the  experi¬ 
mental  results  are  given  in  Section  5.3.  In  Section  5.2,  ijsaediaf "ly 
following,  the  theoretical  basis  for  the  experiments  is  given. 


FIG.  5.4  CRACKED  TEST  SECTION 
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5.2  Theory 


5.2.1  Shield  Velocity  for  Internal  Loading 

For  comparison  with  the  experiments  we  calculate  the  shield 
velocity  for  internal  loading  using  one-dimensional  linear  analysis. 

Figure  5.5  shows  the  analysis  for  a  uniform  load,  taker  so  for  simplicity. 
Upon  arrival  at  the  interface,  the  incident  pulse  in  the  mandrel  (state  1) 
results  in  a  transmitted  shock  in  the  shield  and  a  reflected  rarefaction 
in  the  mandrel,  corresponding  to  state  2.  The  shock  in  the  shield  reflects 
from  the  free  surface  as  a  rarefaction  producing  state  02.  The  arrival 
of  the  rarefaction  at  the  interface  results  in  separation.  Thus,  for 
internal  loading,  the  maximum  duration  of  the  pulse  acting  on  shield  is 
twice  the  transit  time  through  the  shield. 


After  separation,  the  pulse  continues  to  accelerate  the 

mandrel.  As  indicated  in  Fig.  5.5,  if  the  acceleration  continues  long 

enough,  the  mandrel  velocity  (e.g.,  u  _)  eventually  exceeds  the  sepa- 

x  10 

ration  velocity  uQ2  of  the  shield.  To  determine  the  probability  of 

postseparation  impact,  the  motion  of  the  mandrel  and  shield  are  analyzed 
in  Appendix  A.  The  results  show  that,  for  the  thick-walled  mandrel  used, 
postseparation  Impact  does  not  occur. 


Since  separation  occurs  when  the  rarefaction  from  the  free 
surface  arrives  at  the  interface,  the  momentum  1^  transmitted  to  the 
shield  prior  to  separation  is  related  to  the  interface  impulse  I ,  for 
a  pulse  of  arbitrary  shape,  by 


Vi 


(5.1) 


where 


2z 

1  +.  z 
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(a)  x,  t  WAVEFRONT  DIAGRAM 


(b)  p,  u  STATE  DIAGRAM 


FIG.  5.5  x,  t  AND  p,  u  DIAGRAMS  FOR  A  UNIFORM  INTERNAL  PULSE 


and  z  is  the  ratio  of  the  acoustic  impedance*  of  the  shield  to  that 
of  the  mandrel. t 

For  an  incident  pulse  of  arbitrary  shape,  the  mOBenttm  is 

given  by 


I  =  PTv  (5.2) 

where  P  and  T  are  a  characteristic  pressure  and  time  depending  on 
the  pulse  shape,  and  y  is  a  number  depending  on  the  separation  time. 
The  momentum  of  the  shield  is 


1^  ■  phv 


(5.3) 


where  p  and  h  are  the  density  and  thickness  of  the  shield  and  v  is 


* 

Product  of  density  and  sound  speed  (pc). 

For  an  external  load  to  develop  an  interface  pulse  of  duration  com¬ 
parable  to  that  of  the  internal  load  (i.e.,  cot  exceeding  twice  the 
transit  time  through  the  shield),  the  external  load  must  not  exceed  a 
duration  of  twice  the  transit  tloe  through  the  shield.  For  an  external 
pulse  of  arbitrary  shape  satisfying  this  condition,  the  rebound  mo¬ 
mentum  I  /  is  related  to  the  incident  external  impulse  1  by 


I 


I 

e 


where 


Equating  I ' 
pulse 


,  -1 

K  '  * 

r  1  +  z 

to  1^  fro*  (5.1)  yields  for  the  equivalent  external 

1  =  -2T-  I 

e  z  -  1  i 


This  equation  places  no  restrictions  on  pulse  shape,  only  on  impulse. 
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the  velocity.  Putting  (5,3)  and  (5.2)  into  (5.1)  gives 


v 


K. 


? _ T_ 

pc  h/c 


(5.4) 


c  being  sound  speed  and  h/c  the  transit  time  through  the  shield. 
Furthermore,  for  equal  shield  transit  times  h/c,  the  separation  times 
and  hence  the  values  of  y  are  equal.  In  the  experiments  described 
below,  the  transit  tine  of  one  of  the  aluminum  rings  was  equal  to  that 
of  the  phenolic  ring.  For  this  condition  (5.4)  yields 


‘V» 


(pc)* 


(5.5) 


where  v  is  the  observed  velocity  of  the  1/4-inch  aluminum  segment 
and  v  is  the  predicted  velocity  of  the  phenolic  segment.  Values  of 

N 

v  from  (3.3)  are  compared  with  observed  phenolic  segment  velocities 
M 

in  the  next  section. 


3. 3  Experiments 

5.3.1  Experimental  Setup 

Figure  5.6  shows  the  steel  mandrel  and  the  external  parts. 
A  new  setup  was  used  for  each  shot.  The  top  and  bottom  1/8-f ncb-thick 
solid  phenolic  section*  served  as  spacers  to  reduce  end  effects.  Crack¬ 
ing  of  the  4-inch  phenolic  test  section  ess  determined  by  posttest 
examination.  The  material  properties  of  the  phenolic  and  aluminum  and 
of  the  steel  mandrel  are  ^iven  in  Table  5.1. 


The  contact  surfaces  of  the  mandrel  were  ground  to  insure 
positive  contact.  The  three  phenolic  rings  sere  installed  using  oil  be¬ 
tween  the  surfaces.  The  segmented  rings  were  attached  with  a  thin  layer 
* 

of  adhesive.  The  small  test  segments  (Fig.  5.2),  used  to  determine 
the  rebound  velocities  of  the  phenolic  and  aluminum,  were  positioned  at 
right  angles  to  the  focal  plane  of  the  camera. 


* 


Eastman  910. 
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DETONATION 
DISK  \ 


r* 


SCXJO 
PHENOLIC 
SPACER 
0.125“  THICK 


SOLID 
PHENOLIC 
TEST 
SECTION 
0.125“  THICK 


ALUMINUM 

SEGMENTED 

RING 

0.375  THICK 


ALUMINUM 

SEGMENTED 

RING 

0.230“  THCK 

PHENOLIC 

SEGMENTED 

RING. 

0.125  THICK 


SOLID 
PHENOLIC 
SPACER 
0.125"  THICK 


VACUUM-  r' 
CHAMBER  i 
BRACKET  j_ 


Table  5.1 


MATERIAL  PROPERTIES 


7~~ 

Material 

Density 

(g/cm3) 

Sound 
Speed 
(cm/F'Sec ) 

Thickness 

(inch) 

Transit 

Time 

(M-sec) 

Phenolic 

1.39 

0.26 

0.125 

1.23 

Aluminum-1  (6061-T6) 

2.77 

0.  51 

0.250 

1.24 

Aluminum-2  (6061-T6) 

2.77 

0.51 

0.375 

1.86 

Steel 

7.87 

0.  50 

1.0 

5.08 

* 

A  12-lnch  length  of  plastic  explosive  was  placed  in 
contact  with  the  mandrel  core  (Fig.  5.6)  and  a  thin  explosive  detonating 
disk  inserted  at  the  top.  The  di sk, detonated  at  the  center  as  shown 
in  Fig.  5.7,  produced  a  radial  detonation  front  which  simultaneously 
initiated  the  main  charge  around  the  circumference. 

♦ 

EL-506D  ,  DuPont. 


FIG.  5.7  DETONATOR  ARRANGEMENT 
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preliminary  calculations  indicated  that  the  aerodynamic 
pressure  acting  on  the  test  segments  after  separation  would  decelerate 
the  sections  by  15  percent  during  the  1-millisecond  coverage  of  the 
framing  camera.  Consequently,  a  vacuum  chamber  was  designed  to  produce 
a  50-mlcron  environment  (Fig.  5.8).  Figure  5.9  is  a  view  of  the  experi¬ 
mental  setup  wit h  three  light  boxes  installed  on  the  front  of  the  chamber. 


FIG.  5.8  VACUUM  CHAMBER 
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FIG.  5.9  EXPERIMENTAL  SETUP 


5.3.2  Shield  Velocities 


The  segment  velocities  are  given  in  Table  5.2.  The 
observed  values  are  those  obtained  from  photographic  records.  The 
calculated  phenolic  velocities  were  obtained  using  the  velocity  of  the 
1/4-inch  aluminum  segments  and  the  properties  of  Table  5.1  in  (5.5). 

The  last  column  gives  the  ratio  of  the  observed  and  calculated  Micarta 
velocities.  The  agreement  is  vlthln  experimental  error  and  is  better 
than  expected.  The  ratio  is  plotted  in  Pig.  5.10. 

5.3.3  Cracking  Velocity  of  the  Phenolic 

Of  the  shots  listed  in  Table  5.2,  only  the  last  resulted 
in  cracking  of  the  phenolic  test  section  (see  Fig.  5.4).  Thus,  the 
cracking  velocity  of  the  test  section  was  between  4.44  and  4.57  cm/msec. 
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Table  5.2 


SHIELD  VELOCITIES 


Shot 

No.* 

Velocities  (cm/msec) 

V  t 

obs 

Observed 

Calculated 

v 

calc 

3/3"  Al. 

1/4"  Al. 

Phenolic 

Phenolic 

1 

(12276) ^ 

1.40 

1.76 

2.30 

2.18 

|Q 

2 

(12226) 

1.86 

2.27 

2.72 

2.82 

EES 

3 

(11970) 

2.37 

2.81 

3.54 

3.48 

1.02 

4 

(12227) 

3.02 

3.26 

3.93 

4.04 

0.97 

5 

(12277) 

3.17 

3.46 

4.44 

4.29 

1.03 

6 

(12278) 

3.48 

3.71 

4.57 

4.60 

0.99 

* 

SRI  ahot  number  noted  in  parentheses. 

Ratio  of  observed  and  calculated  phenolic  velocities. 


FIG.  5.10  PHENOLIC  VELOCITY 
(ralculafed  vs.  observed) 


y 
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From  (2.29),  the  relationship  between  the  cracking  velocity 

v  and  circumferential  cracking  stress  n  for  a  ring  of  a  linear 
cr  ^  0cr 

material  in  uniform  radial  motion  is  given  by 


cr 


8cr 

pc 


where  we  have  substituted  c  /c 

0 

Using  the  properties  of  Table  5.1 


jr- 

~T 


and 
or  phenolic 


(5.6) 

y  is  Poisson's  ratio, 
in  (5.6)  yields 


v  =  4.29  cm/msec 

cr 

This  is  less  than  10  percent  below  the  velocity  of  4.57  cm/msec  which 
produced  cracking.  The  difference  is  well  within  uncertainties  In 
material  properties  and  experimental  error. 

5.3.4  Pulse  Shape 

To  investigate  the  pulse  shape  at  the  interface,  the 
momentum  of  the  two  aluminum  segments  was  equated  to  the  impulse  de¬ 
veloped  at  the  interface  to  obtain 

t 

mom  =  phv  =  J  p(t)dt  (5.7) 

o 

where  p(t)  is  the  interface  pressure  time  profile  being  investigated 
and  t  is  the  separation  time. 

For  a  uniform  pulse,  (5.7)  yields 


phv  =  P  t  (5.8) 

o  s 

where  Pq  is  peak  pressure.  For  a  pulse  with  a  linear  decay,  (5.7) 

gives  .....  ......  1 ■■  ■, 

.1  ">  .  i .  '•  H  ■ 

phv  =  PQtg  j\  -  ^  (tri» Aguiar) 


(5.9) 


where  Pq  is  peak  pressure  and  T  is  the  pulse  duration.  For  a  pulse 
with  an  exponential  decay 


p(t)  =  PQe 


-Xt 


(5.10) 


where  \  is  the  decay  constant,  (5.7)  yields 


phv 


(exponential) 


(5.11) 


The  two  aluminum  segments  of  different  thicknes -es, 
aluminum-1  and  aluminttm-2, .provide  two  sets  of  momentum  and  t  values. 
For  the  rectangular  pulse  (5.8)  each  est  of  values  defines  a  peak  pres¬ 
sure  ,  P  ,  and  P  The  pressures  obtained  using  the  velocities  from 
0 1  0*2  r 

Table  5.2  are  listed  in  Table  5.3  and  sure  shown  in  Fig.  5.11  for  Shot  3. 

The  ratio  of  the  values  range  from  0.80  to  0.94  indicating  that  the  pulse 

* 

decays  roughly  10  to  20  percent  in  1.24  M-sec.  This  time  is  half  the 
separation  time  for  the  phenolic  and  thin  aluminum  segments;  hence  the 
pulse  probably  decays  roughly  20  to  40  percent  from  0  to  2.5  M-sec.  This 
indicates  that  a  decaying  pulse  is  a  better  representation  of  the  pulse 
actually  realized  at  the  interface. 


As  indicated  in  (5.9)  and  (5.11),  two  parameters,  peak 
pressure  and  decay  time,  are  required  to  specify  a  decaying  pulse. 
Using  the  two  sets  of  values  of  momentum  and  separation  time  from  the 
aluminum-1  and  aluminum-2  segments,  the  peak  pressure  and  decay  tia.c 
may  be  found.  For  the  triangular  pulse 


T 


1 

2 


\\ 
h2  v2 


h2  V2 


1 


2 


(5.12) 


* 


Difference  in  tvaypit  times  of  the  aluminum  segments  (see  Table  5.1). 
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FIG.  5.1V  PULSE  SHAPE  DIAGRAM'  FOR  SHOT  3 


and;  using  T  in  (5,9)  ,  the  peak  pressure  may  be  determined.  For  the 
exponential  pulse 

-Xt, 


h2(l  -  e  *) 


h^l  -  e  *) 


(5,13) 


where  \  also  becomes  a  function  of  v  /v  and  must  be  determined  by 

X  « 

trial  and  error.  The  peak  pressure  is  found  from  (5,11), 

The  velocities  of  Table  5.2  are  plotted  in  Fig.  5.12  as 

a  function  of  incident  impulse.  The  pair  of  lines  indicate  a  constant 

velocity  ratio,  v_/v  ,  of  0,87  which  is  the  rms  value  of  the  expert- 

a  x 

mental  values.  The  variation  of  the  experimentally  determined  values 
from  the  mean  ratio  is  attributed  to  experimental  error.  Since,  for 
the  chosen  pulse  shapes,  the  decay  constants  depend  only  on  the  velocity 
ratio,  which  is  approximately  constant,  the  pulse  duration  at  the  seg¬ 
ment-mandrel  interface  must  be  independent  of  the  incident  impulse. 


Hence,  only  peak  pressure  varies  with  the  explosive  thickness.  The 
values  obtained  for  these  conditions  are  given  in  Table  5.3.  The  corres¬ 
ponding  pulses  for  Shot  3  are  plotted  in  Fig.  5.11. 


Figure  5.13  gives  the  calculated  total  impulse  for  the 
exponential  and  linear  decay  curves  based  upon  the  decay  times  and  peak 
pressures  of  Table  5.3.  The  values  increase  with  incident  impulse  as 
shown.  Although  the  impulse  realized  at  the  interface  should  be  at  least 
a  factor  of  1/R  lower  than  the  Incident  impulse  due  to  geometric  diver¬ 
gence,  the  impulse  levels  for  both  the  pulse  profiles  analysed  are 


FIG.  5.13  CALCULATED  INTERFACE  IMPULSE 
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greater.  The  triangular  representation  lies  16  percea shove  jv! 
the  exponential  representation  lies  40  percent  above.  This 
that 

(1)  of  the  two  pulse  profiles  considered,  the  1 i r 
decay  best  represents  the  actual  pulse,  and 

(2)  the  actual  pulse  probably  decays  more  rapidly 
then  the  linear  decay  rodel. 


-re-,  s 

►  t*  f  * 

!©&  r 
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intentionally  left  bleak. 


APPENDIX  A 


ANALYSIS  OF  POSTSEPARATION  MOTION  OF  MANDREL  AND  SHIEID 

✓ 

A. 1  Introduction 

The  calculated  velocities  of  the  phenolic  segments  given  in 
Section  5  are  valid  only  if  there  is  no  postseparation  impact  during 
the  initial  outward  motion.  The  analysis  given  here  shows  that,  for 
the  chosen  conditions  of  the  experiments,  post separation  impact  does 
not  occur  during  the  initial  outward  motion. 

In  the  analysis,  the  phenolic  ring  is  treated  as  a  linear  oscil¬ 
lator,  taking  the  velocity  at  separation  as  the  initial  velocity.  The 
motion  of  the  outer  surface  of  the  mandrel  is  calculated  by  regarding 
a  thin  outer  layer  of  the  mandrel  as  a  linear  oscillator  being  loaded 
by  pulses  determined  by  the  stress  wave  pattern. 

A. 2  Stress  Wave  Calculations 

Figure  A.l  shows  the  x,t  and  u,p  diagrams  for  a  one-dimensional 
laminate  corresponding  to  the  mandrel  and  the  shield.  Since  wave 
velocities  are  not  very  sensitive  to  pressure,  the  wave  pattern  in  cylin¬ 
drical  geometry  would  be  essentially  the  same.  We  regard  the  load  for 
the  one-dimensional  calculation  as  an  equivalent  load  which  produces  the 
same  pulse  at  the  interface  as  would  occur  in  cylindrical  geometry.  The 
motion  of  the  shield  and  the  outer  surface  of  the  mandrel  are  both  taker, 
to  be  proportional  to  the  load  amplitude  at  the  interface.  Thus,  in  cal¬ 
culating  the  motion  of  the  mandrel,  we  depend  on  the  stress  wave  cal¬ 
culation  only  for  the  frequency  at  which  the  pulses  arrive  at  the  inter¬ 
face. 
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(b)  p,  u  STATE  DIAGRAM 
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FIG.  A.1  x,i  AND  p,u  DIAGRAMS  FOR  A  UNIFORM  INTERNAL  PULSE 


83 


In  Fig.  A. 1  the  shield  separates  iron  the  aandrel  at  the  velocity 
u^.  This  is  taken  as  the  initial  velocity  of  the  phenolic  ring.  The 
initial  velocity  of  the  outer  surface  -‘f  the  n&udrel  is  u^.  The  outer 
layer  of  the  t&cdrel,  indicated  by  the  dashed  line  in  Fig.  A. la,  will 
be  accelerated  further  by  pulses  of  amplitude  p_,  p_,  etc.  JAoreover, 
the  dura tie 3  of  the  pulse  is  proportional  to  the  thickness  of  the  outer 
layer,  hence  the  acceleration  13  Independent  of  the  thickness  chosen. 

In  the  calculations  below,  an  outer  layer  of  convenient  thickness  is 
used  to  calculate  the  action  of  the  cuter  surface  of  the  eandroi  under 

the  action  of  tne  successive  pulses  and  the  restraining  hoop  stresses. 

/ 

A. 3  Motion  of  Phenolic  Ring  and  Mandrel  Surface 


The  station  of  the  phenolic  ring  18  nicply 

02  , 

x  «  -  sin  qt 

Q 


(A.  1) 


where  u  is  the  velocity  of  the  phenolic  at  separation  and  q  is 
given  by 


or 


E  and  p  being  Young’s  isodulua  and  density  for  the  phenolic  shield 
and  r  the  »*an  radius  of  the  shield. 


The  aotion  of  the  outer  surface  of  the  eandrcl  Is  given  by 


x  °  —  sin  *  tiO  _  coa  t J  (A. 3) 

b  % 

where  is  the  velocity  of  the  candrel  surface  at  separation,  F(t) 

is  the  forcing  function,  and 


(A. 4) 


E^  and  being  Young’s  ssodulus  and  density  for  the  steel  candrol  and 

the  tseaa  radius  of  the  outer  layer. 


59 


The  duration  of  the  pulses  at  the  outer  layer  is  twice  the  transit  time 
through  the  layer. 

A. 4  Results 

The  notion  of  the  phenolic  ring  and  the  outer  surface  of  the 
mandrel  are  shown  in  Fig.  A. 2  for  a  mandrel  2b. 4  mm  thick.  The  paths 
indicate  that  impact  does  not  occur  during  tne  initial  outwaru  motion. 
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DISTANCE  —  mra 


FIG.  A.2  DISPLACEMENT-TIME  DIAGRAM  OF  PHENOLIC  SHIELD 
AND  THICK  STEEL  MANDREL 


Fis^re  A. 3  shows  the  sane  phenolic  notion  and  tna>v.  of  a  manovel  8  ran 
thick.  In  this  case  the  mandrel  overtakes  the  shie.«.^,. 

In  the  experiments  a  mandrel  thickness  of  25.4  mm  was  used.  Hence, 
po3t separation  impact  during  the  initial  outward  motion  was  avoided. 
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PRESSURE  — ktwr 


FIG.  A.3  DISPLACEMENT-! IMF  DIAGRAM  OF  PHENOLIC  SHIELD 
AND  th:n  STEEL  MANDREL 


»3 


PflESS^»t 


FIG.  B.3  EXPERIMENTAL  SETUP 


Table  B.l 

SLAB  EXPERIMENT  RESULTS 


Shot 

Impulse 

(taps) 

Loading 

Damage 

1 

1 

10,000 

15  rail 
sheet 

solid  HE 

Severe  spall 

2 

1 

5.000 

15  mil 

1/8 -inch 

Moderate  spall 

1/8  inch 

3 

3,400 

15  mil  1/8-incn 
strips  spaced 

1/4  inch 

Marginal  end  crack 

In  the  final  test,  a  3400-tnp  load  produced  a  minor  crack  at  the 
opposite  eud  from  that  at  which  the  explosive  was  initiated  (Fig.  B.6). 
This  damage  is  attributed  to  edge  effects.  The  3400-tap  loading  is 
considered  a  lower  bound  on  the  impulse  to  cause  separation. 
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FIG.  B.4  COMPOSITE  (10,000  Tcr,s) 
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